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' Abstract 

, In this article, we study the internal stabilization and control of the critical nonlinear Klein- 

Gordon equation on 3-D compact manifolds. Under a geometric assumption slightly stronger than 
the classical geometric control condition, we prove exponential decay for some solutions bounded 
p I ' in the energy space but small in a lower norm. The proof combines profile decomposition and 

I microlocal arguments. This profile decomposition, analogous to the one of Bahouri-Gcrard [2] on 

M'^, is performed by taking care of possible geometric effects. It uses some results of S. Ibrahim 
[2T] on the behavior of concentrating waves on manifolds. 
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Introduction 

In this article, we study the internal stabilization and exact controllability for the defocusing critical 
nonlinear Klein-Gordon equation on some compact manifolds. 

, , J Dii = (9|m — An = —u—\u\^u on [0, -|-oo[xM 

^ ^ I {u{0),dtu{0)) = {uo,ui)e£. 

where A is the Laplace-Beltrami operator on M and £ is the energy space H^{M) x L'^(M). The 
solution displays a conserved energy 

(2) ^(^^ = Uf \9M'+[ 1^1'+/ + H'- 

^ \Jm j m Jai J JM 

This problem was already treated in the subcritical case by B. Dehman, G. Lebeau and E. Zuazua 
[11| . The problem is posed in a different geometry but their proof could easily be transposed in our 
setting. Yet, their result fails to apply to the critical problem for two main reasons, as explained in 
their paper : 

(a) The boot-strap argument they employ to improve the regularity of solutions vanishing in the zone 
of control cj so that the existing results on unique continuation apply, does not work for this critical 
exponent. 

(b) They can not use the linearizability results by P. Gerard |19j to deduce that the microlocal defect 
measure for the nonlinear problem propagates as in the linear case. 

In this paper, we propose a strategy to avoid the second difficulty at the cost of an additional 
condition for the subset uj. It was already performed by B. Dehman and P. Gerard |8] in the case of 

with a flat metric. In fact, in that case, this defect of linearisability is described by the profile 
decomposition of H. Baliouri and P. Gerard [2j. The purpose of this paper is to extend a part of this 
proof to the case of a manifold with a variable metric. This more complicated geometry leads to extra 
difficulties, in the profile decomposition and the stabilization argument. We also mention the recent 
result of L. Aloui, S. Ibrahim and K. Nakanishi [1] for . Their method of proof is very different 
and uses Morawetz-type estimates. They obtain uniform exponential decay for a damping around 
spatial infinity for any nonlinearity, provided the solution exists globally. This result is stronger than 
ours, but their method does not seem to apply to the more complicated geometries we deal with. 

We will need some geometrical condition to prove controllability. The first one is the classical 
geometric control condition of Ranch and Taylor |32j and Bardos Lebeau Ranch [3j, while the second 
one is more restrictive. 
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Assumption 0.1 (Geometric Control Condition). There exists Tq > such that every geodesic 
travelling at speed 1 meets uj in a time t < Tq. 

Definition 0.1. We say that {xi,X2,t) G x R is a couple of focus at distance t if the set 

Fxi,x2,t := e SI^M \ expxM = X2] 

of directions of geodesies stemming from xi and reaching X2 in a time t, has a positive surface 
measure. 

We denote T focus the infimum of the t G M such that there exists a couple of focus at distance t. 

If M is compact, we have necessarily Tjqcus > 0. 

Assumption 0.2 (Geometric control before refocusing). The open set oj satisfies the Geometric 
Control Condition in a time Tq < Tf^cus- 

For example, for T^, there is no refocusing and the geometric assumption is the classical Geometric 
Control Condition. Yet, for the sphere S"^, our assumption is stronger. For example, it is fulfilled if 
a; is a neighborhood of {X4 = 0}. We can imagine some geometric situations where the Geometric 
Control Condition is fulfilled while our condition is not, for example if we take only a a neighborhood 
of {x4 = 0, X3 > 0} (see Remark 10. II and Figure [1] for S"^). We do not know if the exponential decay 
is true in this case. 

The main result of this article is the following theorem. 

Theorem 0.1. Let Rq > and uj satisfying Assumption W.SX Then, there exist T > and 6 > 
such that for any {uo,ui) and {uo,ui) in x L^, with 

\\iuo,ui)\\Hi^L2 < Rq; \\{uo,ui)\\hixl^ < ^0 

||('"0,'"l)|lL2xi^-i < 5; \\iuo,Ui)\\L2^H-i < 6 

there exists g G L°°([0, T], L^) supported in [0,T] x uj such that the unique strong solution of 

j n\u + u+\u\^u = g on [0,T] x M 
\ {u{0),dtu{0)) = (no,ni). 

satisfies {u(T) , dtu{T)) = (no,ni). 

Let us discuss the assumptions on the size. In some sense, our theorem is a high frequency 
controllability result and expresses in a rough physical way that we can control some "small noisy 
data". In the subcritical case, two similar kind of results were proved : in Dehman Lebeau Zuazua 
similar results were proved for the nonlinear wave equation but without the smallness assumption in 

X while in Dehman Lebeau |10J, they obtain similar high frequency controllability results for 
the subcritical equation but in a uniform time which is actually the time of linear controllability (see 
also the work of the author |29| for the Schrodinger equation). Actually, this smallness assumption 
is made necessary in our proof because we are not able to prove the following unique continuation 
result. 

Missing theorem. u = is the unique strong solution in the energy space of 

( Du + u + \u\^u = on [0,T]xM 
\ dtu = on [0, T] X UJ. 
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In the subcritical case, this kind of theorem can be proved with Carleman estimates under some 
additional geometrical conditions and once the solution is known to be smooth. Yet, in the critical 
case, we are not able to prove this propagation of regularity. Note also that H. Koch and D. Tataru 
|26| managed to prove some unique continuation result in the critical case, but in the case u = 
on Lo instead of dtu = 0. In the case of with flat metric and lo the complementary of a ball, B. 
Dehman and P. Gerard [8j prove this theorem using the existence of the scattering operator proved 
by K. Nakanishi |31| . which is not available on a manifold. 

Moreover, as in the subcritical case, we do not know if the time of controllability does depend on 
the size of the data. This is actually still an open problem for several nonlinear evolution equations 
such as nonlinear wave or Schrodinger equation (even in the subcritical case). Note that for certain 
nonlinear parabolic equations, it has been proved that we can not have controllability in arbitrary 
short time, see |15) or |14| . 

The strategy for proving Theorem 10.11 consists in proving a stabilization result for a damped 
nonlinear Klein-Gordon equation and then, by a perturbative argument using the linear control, to 
bring the solution to zero once the energy of the solution is small enough. Namely, we prove 

Theorem 0.2. Let Rq > 0, uj satisfying Assumption \0.S\ and a £ C°°{M) satisfying a{x) > r/ > 
for all X £ UJ. Then, there exist C, 7 > and 6 > such that for any {uo,ui) in x L/^ , with 

the unique strong solution of 

, . f Uu + u + \u\'^u + a{xfdtu = on [0, T] x M 

^' I (7x(0),atn(0)) = (nc^xi). 

satisfies E{u){t) < Ce-'''^E{u){0). 

This theorem is false for the classical nonlinear wave equation (see subsection I3.1.ip and it is why 
we have chosen the Klein-Gordon equation instead. 

Let us now discuss the proof of Theorem 10. 2| following B. Dehman and P. Gerard [8] for the case 
of . We have the energy decay 

E{u){t) = E{u){0) - [ [ \a{x)dtu\^. 

Jo Jm 

So, the exponential decay is equivalent to an observability estimate for the nonlinear damped equa- 
tion. We prove it by contradiction. We are led to proving the strong convergence to zero of a 
normalized sequence Un of solutions contradicting observability. In the subcritical case, the argu- 
ment consisted in two steps 

• to prove that the limit is zero by a unique continuation argument 

• to prove that the convergence is actually strong by linearization and linear propagation of 
compactness thanks to microlocal defect measures of P. Gerard [18] and L. Tartar |35j. 

By linearization, we mean (according to the terminology of P. Gerard [19j) that we have 

WWn — Vn\\\ — ^ where f„ is solution of the linear Klein-Gordon equation with same initial data : 

n— >oo 

r Ovn + Vn = on [0, T]x M 

{ {Vn{0),dtVnm = {Un{0),dtUnm. 
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In our case, the smallness assumption in the lower regularity x makes that the limit is 
automatically zero, which allows to skip the first step. In the subcritical case, any sequence weakly 
convergent to zero is linearizable. Yet, for critical nonlinearity, there exists nonlinearizable sequences. 
Hopefully, in the case of M^, this defect can be precisely described. It is linked to the non compact 
action of the invariants of the equation : the dilations and translations. More precisely, the work 
of H. Bahouri and P. Gerard [2J states that any bounded sequence Un of solutions to the nonlinear 
critical wave equation can be decomposed into an infinite sum of : the weak limit of u„, a sequence of 
solutions to the free wave equation and an infinite sum of profiles which are translations-dilations of 
fixed nonlinear solutions. This decomposition was used by the authors of [8| to get the expected result 
in M^. Therefore, we are led to make an analog of this profile deomposition for compact manifolds. 
We begin by the definition of the profiles. 

Definition 0.2. Let Xoo G M and {f,g) £ = {H^ x L^){T^^M). Given [{f,g),h,x] £ S^^ x 
(M^ X M)^ such that limn{hn,Xn) = (0, Xqo) We call the associated concentrating data the class of 
equivalence, modulo sequences convergent to in £, of sequence in £ that take the form 

(4) 

in some coordinate patch Um ^ U C containing x^o and for some ^jj G C^{U) such that 
^u{x) = 1 in a neighborhood of Xoo- (Here, we have identified Xn,Xoo with its image in U). 

We will prove later (Lemma II. 3p that this definition does not depend on the coordinate charts 
and on ^[/: two sequences defined by @ in diff^erent coordinate charts are in the same class. In what 
follows, we will often call concentrating data associated to \{fig),h^x\ an arbitrary sequence in this 
class. 

Definition 0.3. Let (tn) a bounded sequence in M converging to too cin-d {fn,gn) cl concentrating data 
associated to [{f,g),h,x\- A damped linear concentrating wave is a sequence Vn solution of 

f Ovn + Vn + a{x)dtVn = on M X M 

1 {Vn{tn)-,dtVn{tn)) = (/n,5n)- 

The associated damped nonlinear concentrating wave is the sequence Un solution of 
, , J Dti^ + Un + a{x)dtUn + = on M X M 

^' I {un{Q),dtUnm = K(0),9tT;„(0)). 

If a = 0, we will only write linear or nonlinear concentrating wave. 

Energy estimates yields that two representants of the same concentrating data have the same 
associated concentrating wave modulo strong convergence in Lf^^{M,£). This is not obvious for the 
nonlinear evolution but will be a consequence of the study of nonlinear concentrating waves. 

It can be easily seen that this kind of nonlinear solutions are not linearizable. Actually, it can be 
shown that this concentration phenomenon is the only obstacle to linearizability. We begin with the 
linear decomposition. 

Theorem 0.3. Let (vn) be a sequence of solutions to the damped Klein-Gordon equation ^ with 
initial data at time t = Q (^m'lpn) bounded in £. Then, up to extraction, there exist a sequence of 
damped linear concentrating waves (p^-'-'), as defined in Definition \ 0.3\. associated to concentrating 
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data [{(p^^\ip'^^'>),h^^\x^^\t^^^], such that for any I G 
(7) 
(8) 
(9) 



VT > 0, lim 



Vn{t,x) = V{t,x) + +wli\t,x), 

(0 ^0 

L°=([-T,T],L6(Af))nL5([-T,T],LiO) l^oo 



(Vn,dtV, 



(i)^ 



+ 



oo, 



where o(l) is uniform for t G [—T,T]. 



The nonlinear flow map follows this decomposition up to an error term in the strong following 
norm 



\u\ 



\u\ 



L°°{I,m{M)) 



+ 



L°°{I,L^{M)) 



+ u 



L'^{I,L'^0{M)) ■ 



Theorem 0.4. Let T < Tf 

ocus/'^- Let Ufi be the sequence of solutions to damped non linear Klein- 
Gordon equation 1^ with initial data, at time 0, [(pn,ipn) bounded in £. Denote (resp v the 
weak limit) the linear damped concentrating waves given by Theorem 10.31 and qn^ the associated 
nonlinear damped concentrating wave (resp u the associated solution of the nonlinear equation with 
{u,dtu)t=o = {v,dtv)t=o). Then, up to extraction, we have 



(10) 
(11) 



X)=U + Y Qn^ (*' ^) + (t, x) + r^^ 



lim 

n— ^oo 



l-T,T] «->oo 







where is given by Theorem W. 

The same theorem remains true if M is the sphere and a 
any assumption on the time T. 



(undamped equation) without 



The more precise result we get for the sphere will not be useful for the proof of our controlla- 
bility result. Yet, we have chosen to give it because it is the only case where we are able to describe 
what happens when some refocusing occurs. 

This profile decomposition has already been proved for the critical wave equation on by H. 
Bahouri and P. Gerard [2j and on the exterior of a convex obstacle by I. Gallagher and P. Gerard [17]. 
The same decomposition has also been performed for the Schrodinger equation by S. Keraani [25\ and 
quite recently for the wave maps by Krieger and Schlag |27) . Note that such decomposition has proved 
to be useful in different contexts : the understanding of the precise behavior near the threshold for 
well-posedness for focusing nonlinear wave see Kenig Merle [21] and Duyckaerts Merle [13]. the study 
of the compactness of Strichartz estimates and maximizers for Strichartz estimates, (see Keraani 
|25|). global existence for wave maps |27]... May be our decomposition on manifolds could be useful 
in one of these contexts. Let us also mention that, this kind of decomposition appears for a long 
time in the context of Palais-Smale sequences for critical elliptic equation and optimal constant for 
Sobolev embedding, but with a finite number of profiles, see Brezis Coron |;4j, the book [12] and the 
references therein ... 

Let us describe quickly the proof of the decomposition. The linear decomposition of Theorem 
is made in two steps : first, we decompose our sequence in a sum of an infinite number of 
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sequences oscillating at different rate /i„ . Then, for each part oscillating at a fixed rate, we extract 
the possible concentration at certain points. We only have to prove that this process produces a rest 
wl^ that gets smaller in the norm L°°L^ at each stage. Once the linear decomposition is established, 
Theorem 10.41 says, roughly speaking, that the nonlinear flow map acts almost linearly on the linear 
decomposition. To establish the nonlinear decomposition we have to prove that each element of the 
decomposition do not interact with the others. For each element of the linear decomposition, we are 
able to describe the nonlinear solution arizing from this element as initial data. The linear rest wl^ is 
small in L°°([— T, T],L^) for / large enough and so the associated nonlinear solution with same initial 
data is very close to the linear one. The behavior of nonlinear concentrating waves is described in 
[21) (see subsection 12.2.11 for a short review). Before the concentration, linear and nonlinear waves 
are very close. For times close to the time of concentration, the nonlinear rescaled solution behaves 
as if the metric was flat and is subject to the scattering of R^. After concentration, the solution is 
close to a linear concentrating wave but with a new profile obtained by the scattering operator on 

We finish this introduction by a discussion on the geometric conditions we imposed to get our 
main theorem. For the linear wave equation, the controllability is known to be equivalent to the so 
called Geometric Control Condition (Assumption 10. ip . This was first proved by Ranch and Taylor 
[32) in the case of a compact manifold and by Bardos Lebeau Ranch [3j for boundary control (see 
Burq Gerard for the necessity [5j). For the nonlinear subcritical problem, the result of [lOj only 
requires the classical Geometric Control Condition. Our assumption is stronger and we can naturally 
wonder if it is really necessary. It is actually strongly linked with the critical behavior and nonlinear 
concentrating waves. Removing this stronger assumption would require a better understanding of 
the scattering operator of the nonlinear equation on (see Remark 10. ip . However, we think that 
the same result could be obtained with the following weaker assumption. 

Assumption 0.3. uj satisfies the Geometric Control Condition. Moreover, for every couple of focus 
{xi,X2,t) at distance t, according to Definition \0.1l each geodesic starting from x\ in direction ^ such 
that expxitS, = X2 meets uj in a time < s < t. 

Finally, we note that our theorem can easily be extended to the case of with a metric flat at 
infinity. In this case, our stabilization term a{x) should fulfill the both assumptions 

• there exist R > and p > such that a{x) > p for |x| > R 

• a{x) > p for X G a; where uj satisfies Assumption 10.21 

The proof would be very similar. The only difference would come from the fact that the domain is 
not compact. So the profile decomposition would require the "compactness at infinity" (see property 
(1.6) of [2]). Moreover, the equirepartition of the energy could not be made only with measures but 
with an explicit computation (see (3.14) of [8j) 

Remark 0.1. In order to know if our stronger Assumptions \0. 4^ or \0.3\ are really necessary compared 
to the classical Ceometric Control Condition, we need to prove that the following scenario can not 
happen. We take the example of with uj a neighborhood of {x^ = 0, X3 > 0}. 

Take some data concentrating on the north pole, with a Fourier transform (on the tangent plane) 
supported around a direction ^q. The nonlinear solution will propagate linearly as long as it does not 
concentrate : at time t it will be supported in an neighborhood of the point x{t) where x{t) follows the 
geodesic stemming from the north pole at time in direction ^o- Then, if^Q is well chosen, it can avoid 
UJ during that time. Yet, at time ir, the solution will concentrate again in the south pole. According 
to the description of S. Ibrahim f21f . in a short time, the solution will be transformed following the 
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nonlinear scattering operator on R^. 5*0, at time vr + e the solution is close to a linear concentrating 
wave but it concentrates with a new profile which is obtained with the nonlinear scattering operator on 
M^. This operator is strongly nonlinear and we do not know whether the new profile will be supported 
in Fourier near a new direction . // it happens, the solution will then be supported near the point 
y{t) where y[t) follows the geodesic stemming from the south pole at time vr in direction ^i. In this 
situation, it will be possible that the trajectory y{t) still avoids oj. If this phenomenon happens several 
times, we would have a sequence that concentrates periodically on the north and south pole but always 
avoiding the region u (which in that case satisfies Geometric Control Condition). 




Figure 1: Possible situation on the sphere 

We are led to the following informal question. If S is the scattering operator on M^, is that possible 
that for some data {f,g) G x supported in Fourier near a direction ^q, the Fourier transform 
of S{f,g) is supported near another direction ^i. In other words, can the nonlinear wave operator 
change the direction of the light? 

The structure of the article is as follows. The first section contains some preliminaries that will 
be used all along the article : the existence theorem for damped nonlinear equation, the description 
of the main properties of concentrating waves and the useful properties of the scales necessary for the 
linear decomposition. The second section contains the proof of the profile decomposition of Theorem 
10.31 and 10.41 It is naturally divided in two steps corresponding to the linear decomposition and the 
nonlinear one. We close this section by some useful consequences of the decomposition. The third 
section contains the proof of the main theorems : the control and stabilization. 

Note that the main argument for the proof of stabilization is contained in the last section [3] : 
in Proposition 13.11 we apply the linearization argument to get rid of the profiles while Theorem 13.11 
contains the proof of the weak observability estimates. We advise the hurried reader to have a first 
glance at these two proofs in order to understand the global argument. 



0.1 Notation 

For / an interval, denote 

infill/ = \W\\l°^{I,H^M)) + ll^t'"llL°°(/,L2(M)) + II'"IIl5{7,L10{M)) • 

Moreover, when we work in local coordinate, we will need the similar norm (except for instead 
of H^) 

lll^lll7xR3 = ll^llL°°(/,i/i{R3)) + ll^t^llL°°(7,L2(IR3)) + II^IIl5(/,L10(R3)) • 
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Note that the if I = M, |||ii|||j^]g3 is invariant by the translation and scaling u i— )• {^-j^, ^ h" )- 
The energy spaces are denoted 

£ = H^{M) X L2(M) 

with the respective norms 

li(/,<?)ii^ = ii/iii2(M) + iiv/iii.(M) + 

ll(/>5)llL = \\'^f\\h{T.,^M) + \\9\\l^T^^M) 

We will denote (•, and (•, ■)^^ the associated scalar product. 

When dealing with solutions of non linear wave equations on M (or on T^^M), "the unique strong 
solution" will mean the unique solution in the Strichartz space L^^^(M, L^^{M)) (or Lf^^(R, L^^(Tx^M))) 
such that {u,dtu) E C{R,£) (or C{R,£^^). 

All along the article, for a point x G M, we will sometime not distinguish x with its image in 
a coordinate patch and will write instead of T^^M. M will always be smooth, compact and 
the number of coordinate charts we use is always assumed to be finite. We also assume that all the 
charts are relatively compact. In all the article, C will denote any constant, possibly depending on 
the manifold M and the damping function a. We will also write < instead of < C for a constant C. 
i?|g^(M) denotes the Besov space on M defined by 



u 



l[0,l[(V-^Af)w. ^^^_^+SUp l[2fc,2fc+i[(\/~^A/j^t 



H^{M) 



We use the same definition for I?! oo(^ ) with Am replaced by A^s which can be expressed using the 
Fourier transform and Littlewood-Paley decomposition. Of course, -B| ^^(M) is linked with i?| ^^(M^) 
by the expression in coordinate charts. This will be precised in Lemma |2.1[ 

From now on, a = a{x) will always denote a smooth real valued function defined on M. 

1 Preliminaries 
1.1 Existence theorem 

The existence of solutions to our equation is proved using two tools : Strichartz and Morawetz 
estimates. Strichartz estimates take the following form. 

Proposition 1.1 (Strichartz and energy estimates). Let T > and {p,q) satisfying 

1 3 1 

- + - = 7T, P>2 
p q 2 

Then, there exists C > such that any solution u of 

nv + v + a{x)dtv = f on [-T, T] x M 
{v{0),dMO)) = {uo,ui). 

satisfies the estimate 

\\{v,dtv)\\^^(^^_T^j,]^£^ + \\v\\lp(^[_t,t],Li{M)) < C!i\\iuo,Ul)\\£ + ||/|Ili([_t,T],L2))- 
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Proof. The case with a = for the wave equation can be found in L.V. Kapitanski |23| . To treat 
the case of damped Klein-Gordon, we only have to absorb the additional terms and get the desired 
estimate for T small enough. We can then reiterate the operation to get the result for large times. □ 



Then, we are going to prove global existence for the equation 
(12) 



\I\u + u + \u\^u = a{x)dtu + g on [-T,T] x M 
{u{0),dtu{0)) = {uo,ui)e£ 

with g G L^{[-T,T],L^{M)) and a G C°°(M). 

The proof is now very classical, see for example |36| for a survey of the subject. The critical 
defocusing nonlinear wave equation on was proved to be globally well posed by Shatah and 
Struwe [34i i33| | using Morawetz estimates. Later, S. Ibrahim an M. Majdoub managed to apply this 
strategy in the case of variable coefficients in [22J, but without damping and forcing term. In this 
subsection, we extend this strategy to the case with these additional terms. We also refer to the 
appendix of [2] where the computation of Morawetz estimates on is made with a forcing term. 
We also mention the result of N. Burq, G. Lebeau and F. Planchon [7] in the case of 3-D domains. 

We only have to check that the two additional terms do not create any trouble. Actually, the main 
difference is that the energy in the light cones is not decreasing, but it is locally "almost decreasing" 
(see formula (I13p ) and this will be enough to conclude with the same type of arguments. 

As usual in critial problems, the local problem is well understood thanks to Strichartz estimates 
while we have to prove global existence. We only consider Shatah Struwe solutions, that is satisfying 
Strichartz estimates and we wave uniqueness for local solutions in this class. We assume that there is 
a maximal time of existence to and we want to prove that it is infinite . The solution considered will 
be limit of smooth solutions of the nonlinear equation with smoothed initial data and nonlinearity. 
Therefore, the integrations by part are licit by a limiting argument. 

We need some notations. To simplify the notations, the space-time point where we want to extend 
the solution will be zq = {tQ,XQ) = (0,0). ip is the geodesic distance on M to = defined in a 
neighborhood ?7 of 0. Denote for some small a < /3 < 



m 



{z = {t, x) G [a, /3] X U \ f < \t\ } backward truncated cone 
{z = {t, x) £ [a, f3] X U \ip = \t\ } mantle of the truncated cone 

{x G U \ (p < \t\} spacelike section of the cone at time t 



In what follows, the gradient, norm, density are computed with respect to the Riemannian metric on 
M (for example, we have \V(p\ = 1). We also define 

e{u){t,x) := I ^|f?tn|^ + |Vu|^^ + g |m|^ local energy 

E{u, D{t)) := ^D(t) ^(^)(^' x)dx energy at time t in the section of the cone 

Flux{u, M^)) := ^ IdfuS/ip — Vu|'^ + g \u\^ da flux getting out of the truncated cone 

Lemma 1.1. Let u be a solution of equation U^) . The function E{u, D{t)) satisfies for a < P < 

E{u,D{l3))+Flux{u,M^) = E{u,D{a))+ [[ a{x)\dtu\'^ - ^ [ [ udtu + ^jj gdtu 
and it has a left limit in t = as a function of t. 
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Proof. The identity is obtained by multiplying the equation by dtu to get dte{u) — ?R.div{dtuV xu) = 
a{x)\dtu\'^ — udtu + ^gdtu, then, we integrate over the truncated cone Ka and use Stokes formula. 
Denote f{t) = E{u, D{t)). Using the positivity of the flux and Holder inequality, we estimate 

< /(«) + C{fi - a) + C|a|(/3 - a) + \\9\\L\[a,p\,L^) ll/llL-([a,/3]) • 

Using C\a\{l3-a) ([^_^]) < C'(/5-a)(ll/llLi([^,^]) + ll/llL-([a,/3]))' ^e get for /3-a smah enough. 



(13) fiPf'^ < 1 _ _ ^) + - ") + Ml^[W,PW) 

This property will replace the decreasing of the energy that occurs without damping and forcing term 
in all the rest of the proof. It easily implies that / has a left limit. □ 

Lemma 1.2. For u and g a strong solution of 

Uu + \u\'^u = g on [-r,0[xM 

we have the estimate 



1 



a/2 



|n| 

D{a) 



< Ci^ (/(/?)) + /(/3)i/3) + - /(a)| + 11511^1^2(^/3) \\dtn\\^^^,^j,,^ 



a 



1/3 



+ - /(a)l + Ibllili2(^^) l|f^t'«llLooi2(;^^) 



+(/3-a) sup /(t) + /(t)i/3 

te[a,/3] L 

where we have used the notation f[t) = E{u, D{t)). 

Proof. It is a consequence of Morawetz estimates. The only difference is the presence of the forcing 
term g and the metric. The case of flat metric is treated in |2j. The metric leads to the same estimates 
with an additional term (/? — a) sup f{t) + f{t)^^^ as treated in [22j. Another minor difference is 

that in the presence of a forcing term, the energy does not decrease and /(/3) + f{f3)^^'^ have to be 
replaced by the supremum on the interval. Note also that our estimate is made in the backward cone 
while the computation is made in the future cone in these references. We leave the easy modifications 
to the reader. □ 

The previous estimates will be the main tools of the proof. It will be enough to prove some non 
concentration property in the light cone for L°°L^ , L^L^^ and finally in energy space. It is the object 
of the following three corollaries. 

Corollary 1.1. 

\u{a, x)\^dx — > 0. 



/ 
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Proof. We are going to use the previous Lemma [1.21 replacing ghy g — u + a{x)dtu and with /3 = ea, 
< e < 1. Denote L the limit of /(t) as t tends to given by Lemma ll.ll So for a small enough, 
we have for a constant C > 

We also use 

\\g-u + a{x)dtu\\^,^,^j^,^ < \\g\\LiL2^K^) + Ci/3 - a) + C(/3 - a) \\dtu\\^^^,^^0^ 

which tends to as /3 tends to 0. This yields 

Ih^ / \u{a,x)fdx<Ce{L + L^/^). 

JD(a) 

□ 

Corollary 1.2. 

Proof. Localized Strichartz estimates in cones (see Proposition 4.4 of [22]) give 
IMl^l-^^^rO) < C^('u,£'(s))^/^ + Ihlli5iio(^o) + \\a{x)dtu - u + g\\ ^1^2 (^KO) 

< CE{u, D{s))^/^ + \W\\l^l^(^KO) (1 + \Ml4Li2^K0)) + \\dtu\\L^LHKO) + Wah^L^RO) ■ 

A boot-strap argument and Corollary 11.11 give that for s sufficiently close to 0, ||^i||i4j;^i2(;^o) is 
bounded. We get the announced result by interpolation between L'^L^^ and L°°L^. □ 

Corollary 1.3. 

E{u,D{.s))^Q. 

Proof. Let e > 0. Corollarv 11.21 allows to fix s < close to so that ||'u||2,5iio(j^o) < £■ Denote Vs 
the solution to the linear equation 

Uvs + Vs + a{x)dtVs = 0, {vs, dtv)t=s = {u, dtu)t=s 

then, the difference Wg = u — Vg is solution of 

Dws + Ws + a{x)dtWs = -\u\'^u, iws,dtWs)t=s = (0,0). 

Then, for s < t < 0, linear energy estimates give 

Eo{ws,D{t)y/^ < C||n||^5iio(^o) < Ce' 

where we have set 

Eo{w„ D{t)) = [|Vw;,p + \dtws\^] dx. 

2 JD{t,zo) 

Triangular inequality yields 

Eo{u, D{t)f'^ < Eo{v„ D{t)f'^ + Ce\ 

Since Vs is solution of the free damped linear equation, we have Eq{vs, D{t)) — > 0. This yields the 
result with Eq instead of E. The final result is obtained thanks to Corollary 11.11 □ 
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We can now finish the proof of global existence. 

Let e > to be chosen later. By Corollary 11.31 ^i^, D{s)) < e for s close enough to 0. By 
dominated convergence, for any s < close to 0, there exists > so that 

/ e(n)(s) = E{u, D{s, rj)) < 2e 

J ip{x)<to—s+r) 

where E{u, D{s,r])) is the spacelike energy at time s of the cone centered at (to = xq = 0) (see 
Figure [2]). For s close enough to and s < s' < 0, we apply estimate ([13]) in this cone. It gives 

E{u, D{s',v))'/' < C {e{u, D{s, r?))i/2 + |.' - s| + < Ce^^. 

In particular, ||ii||/^oo2,6(^) < Ce^/^ on the truncated cone 

K = {{s',x)\ (p{x) < r] - s' ,s < s' < 0} 
Therefore, choosing e small enough to apply the same proof as Corollary II. 2| we get 

ll«llL5LiO(i^) < +00. 

Since xq = is arbitrary, a compactness argument yields one s < such that II^^II^sq^ o[L1o(M)) ^ +oo. 
Therefore, by Duhamel formula, {u{t), dtu{t)) has a limit in if as t tends to and u can be extended 
for some small t > using local existence theory. 

t ♦ 




Remark 1.1. It is likely that global existence can also be proved using the Kenig-Merle argument 
124^ and the profile decomposition below ( assuming only local existence ) as is done for example in }27}j 
for the wave maps. 

1.2 Concentration waves 

In this section, we give details about concentrating waves that will be useful in the profile decom- 
position. The first lemma states that Definition 10.21 of concentrating data does not depend on the 
choice of coordinate patch and cut-off function ^ij. 

Lemma 1.3. Let [{f,g),h,x\ ^ £ x (M!j_ x Mx)^ such that limn{hn,Xn) = (0, Xqo) then, all the 
sequences defined by formula ^ in different coordinates charts and the cut-off function are 
equivalent, modulo convergence in £. 
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Proof. It is very close to the one of S. Ibrahim |21| where the concentrating data are given in geodesic 
coordinates. So, let Vm ^ V he another coordinate patch and $ : 1/ i— )■ [/ the associated transition 
map. Without loss of generality, we can suppose that Xqo is represented by in [/ and V. We have 
to prove that the sequences 



and 



K'/Hvix)if o D^O), ^g o Z)cD(o)) 



are equivalent in the energy space associated to M or (the volume form and the gradient are not 
the same but the energies are equivalent). By approximation, we can assume (/, 5) G (Cq^(M^))^. We 
make the proof for the part for /, the proof being simpler for g. We remark that the terms coming 
from derivatives hitting on ^>u{x) tend to in L^. Therefore, we have to prove the convergence to 
of 

2 



^u{Hx))D^{x)V f 



K 



^v{x)D^{0)Vf 



D^{0)x - D^{0)Xn 



First, we prove that the cut-off functions ^jj and can be replaced by a unique Let 5 so that 
B{0,6) C V. Let ^' G C^{B{0,S)) such that ^ = 1 in a neighborhood of and has a support 
included in the set of x such that ^y(x) = ^(/(<I>(x)) = 1, so that ^"S/y = ^ and ^{^u o^) = 
Then, on the support of 1 — ^, we have — ^(xn)!! > ^ some e > and some n large enough. 

Therefore, we have 



< CK 



;i - ■^{x))-^u{^{x))D^{x)Vf 



<I>(x) - <I>(x„) 



L2(y) 



V/ 



^(x) - ^(Xn) 



L2(||<I.(x)-<I'(x,0||>e) 



which is for n large enough since / has compact support. Making the same proof for the other 
term, we are led to prove the convergence to of 



(14)< 



h-^ ^{x)D<^{x)Vf 



D^{hnX + Xn)Vf 



^(x) - ^(x„) 
h„. 



^{x)D<S>{0)Vf 



Z)^>(0)x - L>$(0)x^ 

Ti 



^{hnX + Xn) - ^{ 



D$(0)V/(D$(0)x) 



L'2{{x:\x„+hnx\<S} 

By the fundamental theorem of calculus, there exists Zn{x) E [x„, /i„x + x„] such that 
^(hnx+x^)-^(xn) _ |2)(j)(2;„)x| > C \x\ for some uniform C > 0. As V/ is compactly supported, 
we deduce that for |x| large enough, the integral is zero. So, we are led with the norm (jl4p with 
L^(i?(0, C)) instead of i^({x : |x„ + hnx\ < 6}). We conclude by dominated convergence. □ 

Using the previous lemma in geodesic coordinates, we get that our definition of concentrating 
data is the same as Definition 1.2 of S. Ibrahim |21] . 

Remark that for a concentrating data, not be defined invariantly on T^^M, we can 

only define the limit of {xn — Xoo)/hn- The change of coordinates must act on x„ as an element of 
M and not T^^M even if it converges to Xqo- Yet, the functions {f,g) of a concentrating data "live" 
on the tangent space. Moreover, the norm in energy of a concentrating data is the one of its data. 
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Lemma 1.4. Let {un,Vn) a concentrating data associated to [{'f , ip) , h, x} , then, we have 

\\{Un,Vn)\\e = \\{'P,n\s.^+0{l) 

where V^;^ and L'^{Tx^M) are computed with respect to the frozen metric. 

The proof is a direct consequence of Lemma 11.51 and 11.61 below or by a direct computation in 
coordinates. 

The next definition is the tool that will be used to "track" the concentrations. 
Definition 1.1. Let Xoo € M and {f,g) G ^xoo- Given [{f,g),h,x] G S^oo 

X (M:^ X M)^ such that 

limn{hn,Xn) = (0,Xoo)- Let {fn^gn) ^6 0- Sequence hounded in £, we set 

DhMn,gn) {f,g) 

if in some coordinate patch Um ^ U C M.'^ containing x^o CLnd for some '^u G C^{U) such that 
^uix) = 1 in a neighborhood of Xoo, we have 

1 

hn {^ufn, K'^ugn) {xn + Kx) ^ (/, g) Weakly in 

where we have identified ^uifn, g-n) with its representation on T^^M in the local trivialisation. 

If this holds for one ([/, ^'j/), it holds for any other coordinate chart with the induced transition 
map. 

We denote Dj^^fn ^ f if we only consider the first part concerning and Df^^g^ g for the 
L^ part convergence. 

Of course, this definition depends on the core of concentration h and x. In the rest of the paper, 
the rate h and x will always be implicit. When several rate of concentration [hi^\x^^^], j G N, are 
used in a proof, we use the notation D^j^^ to distinguish them. 

The fact that this definition is independent of the choice of a coordinate chart can be seen with 
the following lemma which will also be useful afterward. 

Lemma 1.5. -D/i„(/n)5n) ^ if^g) ^-^ equivalent to 

^Mfn-^MUn > / ^x^f '^x^^ 

gnVn — > I gip 

where {un,Vn) is any concentrating data associated with [{^,il^),h,x\. 

The V is computed with respect to the metric on M when the integral is over M and with respect 
to the frozen metric in Xoo when the integral is over Tx^M. 

Proof. We only compute the first term for the norm and assume ip G C^(M^). duj{y) denotes 
the Riemannian volume form at the point y, -y the scalar product at the point y and 'Vh„x+xri — 
g{hnX + x„)"^V. 
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We denote Vh = -i-/^ and L„y = hi Jy V^^ [^vfn {xn + hnx)] ■ V^^ip{x) duj{0). 



Ln,V = K I Va:„+h„x [^vfn {Xn + Kx)] -(xr^+hnx) ^x„+h„xV{x) duj{Xn + Kx) + o(l) 

= hi ^v{Xn + Kx)iV x„+h„xfn) (Xn + Kx) -{x^+h^x) ^x„+h„x^ix) duj{Xn + h^x) + o(l) 
= hn~' j^Vyfn{y)-y^viy)^y^l^^^^) du{y) + o{l) 



IV 



'^v{y)f 



n 

y-x 

hr. 



du{y) + 0(1) = / VM/n • VA/n„ + o(l) 
JM 



Therefore, Lny tends to J V f{x)-'VLp{x) duj{0) if and only if ^ m fn'^ MUn has the same hmit. □ 

An easy consequence of this lemma is the link with concentrating waves. 
Lemma 1.6. Let {fn,gn) be some concentrating data associated with [{f,g),h,x}, then, we have 

DhMn,9n) {f,g) 

Proof. Lemma 11.31 permits to work in geodesic coordinates so that the metric g is the identity at 

_ 1 

the point Xoo- In this chart, we have /„ (x„ + hnx) = ^u{xn + hnx)hn ^ f ■ So, the computation of 
Lemma [T31 gives / Vqo/ • V oo^dijj{Q) = Jj^jV Mfn • ^MUn + o(l) which gives the result. □ 

We conclude this subsection by a definition of orthogonality that will discriminate concentrating 
data. 

Definition 1.2. We say that two sequences \h!-'^\x^^\t^^^ and \h^'^\ x^'^\ t^'^^ are orthogonal if either 



log 



^1) 



u{2) 



-00 



(1) / (2) 
" •^00 7^ Xoo 



hn^ = hn^ = h and x2d^ = = Xnr, and in some coordinate chart around Xno, we have 



.(1) _ .(2) 



+ 



„(1) ^(2) 



+00 



h h h-^o 

We note \h^^\x(^Kt(^^] ± [h^^\x^^\t(^^ and ±h {x^^Kt^^^) if h^^^ = h^^^ = h. 

This definition does not depend on the coordinate chart. This can be seen because we have the 



estimate ^ 



(1) (2) 



< 



(2). 



< c 



^(1) _ J2) 



if $ is the transition map. 



1.3 Scales 

In this subsection, we precise a few facts that will be useful in the first part of the proof of linear 
profile decomposition which consists of the extration of the scales of oscillation hn. 

On the Hilbert space E = H^{M) x Lp'{M)^ we define the self-adjoint operator Am by : 

D{Am) = Hi X H\, 
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We define similarly A^d with the flat laplacian. We denote A^d the obvious operator on {H^ (M*^) x 
L^(M'^))^ obtained by applying A^d on each "coordinate". 

The following definition is taken from Gallagher-Gerard[T7|. 

Definition 1.3. Let A be a selfadjoint (unbounded) operator on a Hilbert space H. Let (hn) a 
sequence of positive numbers converging to 0. A bounded sequence in H is said {hn)- oscillatory 
with respect to A if 



0. 



(15) lim 
(un) is said stricly [hn)- oscillatory with respect to A if it satisfies M^fi and 



lim 

n— >oo 



0. 



H 

At the contrary, (un) is said {hn)-singular with respect to A if we have 



for all < a < 6. 



Remark that 1 



|x|<l 



can easily be replaced by a well chosen function ip S C^f 



Moreover, if a 



sequence is stricly (/i„)-oscillatory while a second sequence [vn) is (/i„)-singular, then we have 
the interesting property that {un,Vn)fj — > 0. 

n— 5>oo 

Proposition 1.2. Let M = U^-^Ui a finite covering of M with some associated local coordinate patch 
<I>i : C/j — 7- Vi, C M^. Let 1 = ^'j be an associated partition of the unity of M with ^'j G C^{Ui). Let 
{un,Vn) a bounded sequence in the M energy space and hn a sequence converging to 0. Then {un,Vn) 
is (stricly) (hn)- oscillatory with respect to Am, if and only if all the <l>j*^j(u„, are (strictly) 
(hn) -oscillatory with respect to A-^d. 

Proof. First, we remark that a sequence is (strictly) (/i„)-oscillatory with respect to A if and only if 
it is (strictly) (/i^)-oscillatory with respect to A^. So we can replace Am and A-^3 by —{Am,^m) 
and — (A]r3, Aj^a). We apply a proposition taken from |17| that makes the link between oscillation 
with different operators. 

Proposition 1.3 (Proposition 2.2.3 of [17J). Let A : Hi — )■ H2 be a continuous linear map between 
Hilbert spaces Hi, H2- Let Ai be a selfadjoint operator on Hi, A2 be a selfadjoint operator on 
H2. Assume there exists C > such that A{D{Ai)) C D{A2), A*{D{A2)) C D{Ai) and for any 
u G D{Ai), V G D{A2), 



(16) 
(17) 



\\A2Au\\ <Ci\\Aiu\\ + \\u\ 
\AiA*v\\ < C{\\A2v\\ + \\v\\ 



If a bounded sequence [un) in Hi is (strictly) (hn) -oscillatory with respect to Ai, then (Aun) is 
(strictly) (hn) -oscillatory with respect to A2- 

To prove the first implication, we apply the proposition with 
A{u,v) = (<l>i*^i(n, t'), • • • j^Nif'^ n{u,v)). We only prove the necessary estimates, the inclusions of 
domains being a direct consequence of the inequalities and of the density of smooth functions. To 
simplify the notation, we denote {ui,Vi) = ^i.^'^i{u,v). The proof of (jl6p mainly uses the equivalent 
definitions of the i/* norm on a manifold. 



I {Ui, Vi) I 



< 

< 



l^^RSUiW^i + \\Aj^3Vi\\^2 < WuiW^S^ + \\Vi\\H\ 

|n||^i -^\\Amu\\jji +111^11^2 +||Aa/u||^2 

\^Miu,v)\\rrl y r2 + \\(u,v) 



< 



m H-3 + \\v\\fT2 
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Let us prove (|17p for the duality x of the scalar product. Let {f,g) = {fi,gi)i=i—N G 
(C^(M3) X C^(M3))^ and {u,v) G C°^(M). 



{{u,v),AMA*{f,g))Hi(M)xL^M) = (A^m(«>''^)i (/>5))(Hi(R3) 
i i 



xL2(R3))JV 



+ ll(/,^/)ll(/^l3XL^3)- 



and 



Therefore, we get \\AlA* (f, g)\\^,^^^,^ < C f\\Al,^^if,g) 

\ ^ h3 r3'' ^ 

Proposition 11.31 implies that (strict) (/i„)-oscillation of {un) with respect to Am implies (strict) (hn)- 
oscillation of Aun with respect to A^3 ^y. 

To prove the other implication, we use a quite similar operator. Denote ipi some other cut-off 
functions in C^{Vi) C C^(M^) such that c^j = 1 on Supp{^ij,^i). We define T the bounded operator 
from {H^i X L^g)^ to x L^^ given by 

r(/,5) = E^»Vi(/.,50 



Then, we have T o A = Id and we only have to prove that (strict) (/i„)-oscillation of {fn,gn) with re- 
spect to A^3 ]\[ implies (strict) (/i„)-oscillation oiT{fn, gn) with respect to Am- The needed estimates 
are quite similar and we omit them. □ 

Remark 1.2. Another way to prove Proposition {T7^ would have been to use the pseudodifferential 
operators (/?(/i^Am) cls in fB^. 

Now, we will prove that the (/in)-oscillation is conserved by the equation, even with a damping 
term. 

Proposition 1.4. Let T > 0. Let {(pnjipn) o, bounded sequence of £ that is (stricly) (hn) -oscillatory 
with respect to Am- Let Un be the solution of 



(18) 



OUn + Un = a{x)dtUn on [0, T] X M 
{Un{0),dtUniO)) = {ipn,1pn)- 



Then, {un{t) , dtUn{t)) are (strictly) (hn)-oscillatory with respect to Am, uniformly on [0,r]. 

At the contrary, if {fmipn) is (hn)-singular with respect to Am, {un{t) , dtUn{t)) is (hn)-singular with 

respect to Am, uniformly on [0,T]. 

Proof Let x £ ^(r'(^^) such that < x(s) ^ 1 a^nd x(s) = 1 for \s\ < 1. The (/i„)-oscillation (resp 
strict oscillation) is equivalent to lim (1 — x){Fl?h'^I^){un,dtU„ 



(resp lim x(%^)(^in, 5tn„,) 



e R-^oo 



0). 
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(1 — x){R 'T'n^)'^n is solution of 



{Vn{<^).dtVnm 



(19) 

and energy estimates give 



a{x)dtvn - [x{R^hlA),a]dtUn on [0, T] x M 
(l-x)(i?'/i^A)(v9„,V„). 



2 1,2 



lLi([0,t],L2) 



< Ct II (1 - x){R"K^){Vn,i^n)\\s + CtRK. 

where the last inequality comes from the fact that x{~h'^^) is a semiclassical pseudodifferential 
operator, as proved in Burq, Gerard and Tzvetkov [6], Proposition 2.1 using the Helffer-Sjostrand 
formula. 

Therefore, passing to the limitsup in n and using the oscillation assumption, we get the expected 
result uniformly in t for < t < T. The results for strict oscillation and singularity are proved 
similarly. □ 

Proposition 1.5. There exists Ct > such that for every ((/?«, V'n) bounded sequence of £ weakly 
convergent to 0, we have the estimate 

lim^ IIK,5tn„)||^^([0^5.]^5i^^(^^)xB0^(Af)) < ^T^im^ IKv'n, ^n)|lBi_^(M)xBO^{Af) 
where Un is the solution of il8\) . 



Proof. Without loss of generality and since the equation is linear, we can assume that || ((/?„, ^/;„)||£- 
is bounded by 1. Let e > 0. Let xo,X G C^(M) so that 1 = xo + EfcLi x(2"^''x). We denote 



x{2 A)un- Using the same estimates as in the previous lemma, we get 



[u'n{t),dtu'^{t)) 



< Ct 



+ Ct2- 



Take K large enough so that (7^2 ^ < e for k > K so that we have . 



(20) 



{U^{t),dtu'^{t)) < CT\\{^n,'>Pn)\\Bl^(M)xBO^{M) + 



£. 



Then, for k < K, we use again some energy estimates for the equation verified by u^, we get 



{u^nit),dtuUt)) 



< Ct 



+ Ct 



[a,x(-2-''=A)]atu 



Li([0,r],L2) 



Yet, for fixed k, [a,x(— 2 ^'^A)] is an operator from into (for instance) and we conclude by 
the Aubin-Lions Lemma that for fixed k < K 



(21) 



lim 

n— >oo 



(<(t),9t<(t)) 



< Ct lim 



We get the expected result with an additional £ by combining (|2U|) and (|21|) . 



□ 



We end this subsection by two lemma that will be useful in the nonlinear decomposition. The 
first one is lemma 3.2 of 1171. 



Lemma 1.7. Let hn and hn he two orhogonal scales, and let (fn) and fn be two sequences such that 
such Vfn (resp Vfn)) is strictly (hn) (resp hn)-oscillatory with respect to Ajga. Then, we have : 



lim 



fnfn 







19 



Then, we easily deduce the following result. 

Lemma 1.8. Let hn and hn be two orhogonal scales and Vn, Vn he two sequences that are strictly hn 
(resp hn) oscillatory with respect to Am (considered on the Hilbert space H^), uniformly on [—T,T]. 
Then, we have 

\\VnVn\\L--{[-T,TlL--^iM)) 

Moreover, the same result remains true if Vn is a constant sequence v S and hn = 1- 
Proof. Using a partition of unity 1 = ^? adapted to coordinate charts, we have to compute 

||<I>j*^'jU„^'j7}„|| L°°{[-T,T],L'-'{R^)) • 

Using Proposition 11.21 we infer that <I>j*4'jt'„ is strictly (/i„)-oscillatory with respect to A]g3 (defined 
on H^) and the same result holds for V (<^j*^'jt'n) with respect to Aga defined on L^. We conclude 
by applying Lemma 11.71 to ^i^:^iVn and ^i^:^iVn- Q 

1.4 Microlocal defect measure and energy 

In this subsection, we state without proof the propagation of the measure for the damped wave 
equation. We refer to [T8j for the definition and to [19j Section 4 or [16] in the specific context of the 
wave equation). It will be used several times in the article. 

Lemma 1.9. [Measure for the damped equation and equicontinuity of the energy] Let Un, Un be two 
sequences of solution to 

\JUn + Un = a{x)dtUn, 

weakly convergent to in S. Then, there exists a subsequence (still denoted Un, Un) such that for 
any t G [0,T] there exists a (nonnegative if Un = Un) Radon measure ^* on S*M such that for any 
classical pseudodifferential operator B of order 0, we have with a uniform convergence in t 

(22) (i?(-A)i/V(t),(-A)i/2n„(t)) +^BdtUn{t),dtUn{t))LHM) ^ f MB)dfi'. 



Moreover, one can decompose 



fi' = + fit) 



which satisfy the following transport equation 

dtlJ'±{t) = ±H\i:\^i^±{t) + a{x)fi±(t). 

Furthermore, if tn — > t, we have the same convergence with t replaced by tn in h22(l. 

The microlocal defect measure of a concentrating data [{'^,ilj),h,x\ can be explicitely computed, 
as follows 

, r+°° - 2 „ 

/i± = (27r) (x) / ^{ri)±i\ri\^^{ri) r^dr. 

Jo 

This can be easily computed, for instance, with the next lemma. 
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Lemma 1.10. Let {(pn,il^n) = [(V') "0);/^;^] ^6 ^ concentration data and A{x,Dx),B{x,Dx) two 
polyhomogeneous pseudodifferential operators of respective order 0. Then 

\\{A{x, Dx)^Pn, B{x, Dx)'ipn) - [{Ao{xoo, Dx)^p, Bo{xoo, Dx)ip), h, x]\\Hiy^L2 — > 

n— >oo 

where Aq{xoo, D^) is the Fourier multiplier of homogeneous symbol ao{xoo,S,) defined on T*^M 

Proof. We only give a sketch of proof for B{x, Dx)ipn- By approximation, we can assume that 
'ip E Cq°(M^\0). In local coordinates centered at Xqo = 0, we have for a o(l) small in 



B{x,Dx)i>n = hn~'B{x,Dx) 



^u{x){xi') 



+ 0(1) = /i„ ^ Dy)^l^] ( ) + o(l) 



where Bn{y,Dy) is the operator of symbol bn{y,0 — ^oihuV + Xn,C/hn)- Here 6o is the principal 
symbol of homogeneous for large ^. We write 6o(^ny+a;n, ?/^n) = bo{xn,£./hn)+hny /o^(<9y6o)(a^n+ 
thn,£,/hn)dt. The first term converges to 6o(0,^) by homogeneity while the second produces a term 
small in L^. □ 

The previous lemma is made interesting when combined with the propagation of microlocal defect 
measure. 

Lemma 1.11. Let Un a sequence of solutions of Oun + Un = a{x)dtUn weakly convergent to and 
Pn = [('/^) V')iZi)3iii] linear damped concentrating wave. We assume Dfi{un,dtUn) 0. Then, for 
any classical pseudodifferential operators A{x,Dx) of order 0, we have uniformly for t G [—T,T] 

In particular, we have 

Vp„ • Vun + dtpndtun ^ in V'{] - T, T[xM). 
Proof. We first check the property for t = t^- Using Lemma [1.101 several times, we are led to estimate 

where {ipn,ipn) are the concentrating data associated with [{A{xoo, L)x)(p, B{xoo, L)x)ip),h,x}. Then, 
the hypotheses Dh{un,dtUn) and Lemma [1.51 vields the convergence to for this particular case 
t = tn- We conclude by equicontinuity and by the propagation of joint measures stated in Lemma 

\m □ 



2 Profile Decomposition 
2.1 Linear profile decomposition 

The main purpose of this section is to establish Theorem 10.31 It is completed in two main steps : the 
first one is the extraction of the scales /i„ where we decompose w„ in an infinite sum of sequence 
Vn which are respectively /i„ -oscillatory and the second steps consists in decomposing each vii in 
an infinite sum of concentrating wave at the rate hn . Actually, in order to perform the nonlinear 
decomposition, we will need that, in some sense, each profile of the decomposition do not interact 
with the other. It is stated in this orthogonality result. 

Theorem 10. 3f . With the notation of Theorem \U. 'A we have the additional following properties. 
If2T < Tfocus, we have (M'^) , x^^) , t W ) _L {h^^\ x^^\ t^^'^) for any j / k, according to DefinitionlTE 
If M = and a = (undamped solutions), but with T eventually large, we have 
{h!:^\{-l)'^x^^\-^^^ + m-K) orthogonal to {hS^\x^^\^^^) for any m£Z and j / k. 
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2.1.1 Extraction of scales 

Proposition 2.1. Let T > 0. Let [ipn,ipn) ol hounded sequence of £ and Vn the solution of 



(23) 



+ Vn = a{x)dtVn Oil [-T, T] X M 
(u„(0),5tf„(0)) = (v7„,V'n)- 



Then, up to an extraction, Vn can be decomposed in the following way : for any I G N* 

I 



Vn{t,x) = v{t,x) + Y^vi^\t,x) + p^J,\t,x), 



where Vn is a strictly (hn )- oscillatory solution of the damped linear wave equation i2!J]) on M. The 

(i) (?) 
scales hn satisfy hn — > and are orthogonal : 



(24) 



Moreover, we have 



log 



(25) 



lim 

ra— >oo 



L°°(]-T,T[,L6(M)) l- 



l^n 







+00 if j / k. 



(26) ||(^n,5*^;„)(t)||| = ||(z;,a*.;)(t)||| + ^ (^(f-),a*T;(^-))(t) + (p«,9i/>«)(t) +o(l)(t), 
where o(l)(t) — > uniformly for t G [— T, T]. 

n— ^-oo 

Proof. We first make this decomposition for tiie initial data as done in j20j (see also |2j). Then, using 
the propagation of (/in)-oscillation proved in Proposition ll.4[ we extend it for all time. 

More precisely, by applying the same procedure as in |20j, with the operator Am, we decompose 



where {fn\ipn^) is (/in '')-oscillatory for Am, hn' — > 0, and 

n— ^oo 



Si) 



(27) 



lim sup 



l[2^2^+M(^M)($i^*^^): 



E l—^oo 



0. 



Moreover, we have the orthogonality property : 

I 



||(V.n,Vn)|lI = ||(v^,V')llI + E k^n^^-^) + {'^^\^^^) +0{l), n 



J 



oo 



and the hn^ are orthogonal each other as in ([2411 . Moreover, is {h^i^ )- singular for 1 < j < /. 

This decomposition for the initial data can be extended to the solution by 



Vnit,x) = vit,x) +^v';^\t,x) + p'tl\t,x), 
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where each vii is solution of 



a{x)dtv':i^ on x M 
0') 



Thanks to Proposition 11.41 each {vn \t),dtVn' (t)) is strictly (/ijf '')-oscillatory and {pn' {t),dtPn' (t)) 
is (hn^)- singular for 1 < i < /. So, we easily infer for instance that 
{p'n\t),dtPn\t)), {v^n \t),dtvi^\t))) — > Uniformly on \-T, T] where ( , )^ is the scalar product 

on £. This is also true for the product between Vn and Vn , j ^ k thanks to the orthogonality (I24p . 



(i)^ 



The same convergence holds for the product with v by weak convergence to of the other terms. 
Then, we get 



n — >■ oo. 



Let us now prove estimate (|25p of the remaining term in L°°(L®). (|27p gives the convergence to 
zero of (pi'^(O), 5t/3n^(0)) in i?|j^(M) x B2^{M). We extend this convergence for all time with 
Proposition 11.51 and get 



sup lim 



pii\t),dtp'i\t) 



z.oc 2,00 



0. 



The following lemma will transfer this information in local charts. 
Lemma 2.1. There exists C > such that 
1 



c 
1 



^ 11^/11 Bi,^(IR3)iv 



< 



^(M) ^ ||A/||^l^^(]u3)iV 



where A is i/ie operator described in Proposition \l.i^ of cut-off and transition in N local charts. 

We postpone the proof of this lemma and continue the proof of the proposition. Using this lemma, 
we get for every coordinate patch and £ Cg°(C/i). 



lim 

n— >cxD 



L°-{[-T,T],Bl 







The refined Sobolev estimate. Lemma 3.5 of [2], yields for any / E 



L6 



< 



1/3 
L2 



2/3 



< 



1/3 



2/3 

-B^cx) 



Therefore, we have 



and finally 



lim 

n— ^-oo 



lim 

n— >oo 



^(0 







L°°([-T,T],L6(M)) i^oo 
This completes the proof of Proposition 12.11 up to the proof of Lemma 12.1 



□ 
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Proof of Lemma \2. 1\ We essentially use the following fact : see Lemma 3.1 of [2]. Let /„ be a 
sequence of L^(M^) weakly convergent to and compact at infinity 



lim / |/(x)r dx — > 0. 

Then, tends to in i?2oo(^^) ^^^'^ only if is /i„ singular for every scale 

Actually, the same result holds Am, with the same demonstration. The compactness at infinity 
in is only assumed to ensure 



lim 



hoo 



l[-A,A]{A^3)/n 



= for any A > 

L2 



which is obvious in the case of Am because of weak convergence and discret spectrum. 

Using Proposition 11.21 we obtain that /„ is (/i„)-singular with respect to Am if and only if A/„ is 
(/in) -singular with respect to A^3. Combining both previous results, we obtain that the two norms 
we consider have the same converging sequences and are therefore equivalent. □ 

2.1.2 Description of linear concentrating waves (after S. Ibrahim) 

In this subsection, we describe the assymptotic behavior of linear concentrating waves as described 
in [21] of S. Ibrahim. In [2T], it is stated for the linear wave equation without damping. We give 
some sketch of proof when necessary to emphasize the tiny modifications. 

The following lemma yields that for times close to concentration, the linear damped concentrating 
wave is close to the solution of the wave equation with flat metric and without damping. It is Lemma 
2.2 of [21], except that there is an additional damping term which disappears after rescaling. We do 
not give the proof and refer to the more complicated nonlinear case (see estimate (|53l) ). 



Lemma 2.2. Let Vn = [{'■P,'4'),h, x,t)] be a linear damped concentrating wave and v solution of 

OooV = on R X Ta:^M 



(28) 



[V 



Denote Vn the rescaled function associated to v, that is Vn = ^*'^^^^v (^^jf^, ^ h^" ^ where (U,^) is 
a coordinate chart around Xoo and ^ G C^{U) is constant equal to 1 around Xqo- Then, we have 



}^\\\Vn - Vn\\\lt„_AhMAhu]xM ^"^0 

Corollary 2.1. With the notation of the Lemma, if tn = t„ + (C + o(l))/i„, then (^;„, dtVn)\i^i^ is a 
concentrating data associated with [{v{C),dtv),h,x\. 

Moreover, Lemma 2.3 of S. Ibrahim [21j yields the "non reconcentration" property for linear 
concentrating waves. 

Lemma 2.3. Letv = [{(Pjip^hi ^jt)] be a linear (possibly damped) concentrating wave. Consider the 
interval [—T,T] containing too, satisfying the following non-focusing property (see Definition \0. 1\) 

(29) rues (Fx,zoo,s) =0 Va; G M and s / such that t^o + s £ [-T, T]. 

Then, if we set In^ = [—T, tn — A/i„] and In^ =]tn + A/i„, T], we have 

l™||t^n|li^(4,Au/3.A^^6(M)) ^ 
1™ ll^"llL5(,l>Au/3.A^ilO(M)) ^ 0- 
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Sketch of the proof of Le'mma \2.3\ in the damped case. To simplify the notation, we can assume tn = 
0. In [2l] . the proof is made by contradiction, assuming the existence of a subsequence (still denoted 
Vn) such that ||fn('Sn)||j^6(jv/) C > and oo. If s„ — )• r 7^ 0, using Concentration- 

Compactness principle of [30j , we are led to prove that the microlocal defect measure /i associated to 
Vn{tn) satisfies ^({y} x 5^) = for any y G M. We use the same argument for the damped equation 
except that in that case, the measure (J' associated to Vn{t) is not solution of the exact transport 
equation but of a damped transport equation (see for Lemma [1.9p . Yet, the non focusing assumption 
(j29p still implies /i*({y} x S^) = for all y G M and t 7^ 0, which allows to conclude similarly. 

In the case r = 0, we use in local coordinates the rescaled function t^y^(s, if) — ^J~^n^n{,^n^^ ^nU'^'^n) • 
Vn at time s = is a concentrating data at scale hn/ Sn- We prove lim ||{;„(1, •)||^6(jr3) = 0. Again 
by concentration compactness, it is enough to prove that the microlocal defect measure /i* of Vn 
propagates along the curves of the hamiltonian flow with constant coefficient H\^\. Since Vn is 
solution of + Vn + a{x)dtVn = 0, Vn is solution of Dn^n + s^Vn + Sna{sn • +Xn)dtVn = where 
□„ is a suitably rescaled D'Alembert operator. Since the additional terms s^Vn + Sna{sn • +Xn)dtVn 
converges to in L^L^, we can flnish the proof as in Lemma 2.3 of [21] by proving that /i* propagates 
as if □„ was replaced by Doo, that is along the hamiltonian H\^\. 

The estimate in norm L^L^^ is obtained by interpolation of L°°L^ with another bounded Strichartz 
norm. □ 

In the specific case of 5^, Lemma 4.2 of |21) allows to describe precisely the behavior of concen- 
trating wave for large times, as follows. 

Lemma 2.4. Let p he a sequence of solutions of 

( Upn = on [0,+oo[xM 

\ (Pn(0),9tP„(0)) = i(fn,'4^n) 

where (ipnji^n) is weakly convergent to (0,0) in E . Then, we have 

Pn{t + vr, x) = -Pn{t, -X) + o(l)(i) 

where the o(l)(t) is small in the energy space. The same holds for solutions ofOun + Un- 

In particular, ifp is a concentrating wave associated with data [{(p,il^,h,x,t)], then, for any j G N, 
Pn{t + JT^^x) is a linear concentrating wave associated with [{— {^,'<p){{—^y ■),]!,{ — '^y x,t)\ 

In the previous lemma, —x refers to the embedding of into M^. Moreover, the notation 
(93, '(/')(— •) could be written more rigorously {ip,ip)(DcxiT) where DoqI is the differential at the point 
Xoo of the application I : x ^ —x defined from into itself. Actually, we are identifying the tangent 
plane at the south pole with the one on the north pole by the application x 1— )• —x on M^. 

The fact that the result remains true for the equation + u = Q comes from the fact that for 
initial data weakly convergent to zero, the solutions of = and + t; = with same data 
are asymptotically close in the energy space. This can be proved by observing that for a weakly 
convergent sequence of solutions u„ the Aubin-Lions Lemma yields that Un converges strongly to 
in L°°{[—T,T],L?'). So rn = Un — Vn is solution of Dr^ = Un and converges strongly in £. 

2.1.3 Extraction of times and cores of concentration 

In this section, /i„ is a fixed sequence in converging to 0. For simplicity, we will denote it by h 
and Uh for sequences of functions. The main purpose of this subsection is the proof of the following 
proposition, which is the profile decomposition for /i-oscillatory sequences. It easily implies Theorem 
10.31 when combined with Proposition 12.11 
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Proposition 2.2. Let (uh) he a h-oscillatory sequence of solutions to the damped Klein-Gordon 
equation l[23\). Then, up to extraction, there exist damped linear concentrating waves p^, as defined 
in Definition \0.3l associated to concentrating data [{^^''\ip^^'^),h,x^''\t^''^], such that for any I € N* , 
and up to a subsequence, 

(30) Vh{t,x) = ZLlPn\t,x)+W^\t,x), 



(31) VT > 0, lim 



71— >00 



.,(0 



(32) \\ivk,dtv,)re = EU iPh^dtp)i') + {w):',dtw):>) +o(l), as/i^oo 



— > 

L^Q-T,T[,L^{M)) «-s>oo 
(0 «„„(0> ^ 



£ 



£ 

uniformly for t G [—T,T]. 



Moreover, if 2T < Tfocus, for any j ^ k, we have {xj^\tl-''^) _Lft, {x^-'\t^-'^) according to Definition 

TM 

If M = and a = (undamped solutions), but with T eventually large, ((— l)™x*-'^\ t^'^-' + mvr) is 
orthogonal to (x^-'^t *■■'-') for any m € Z and j ^ k. 

Remark 2.1. The assumptions to get the orthogonality of the cores of concentration are related to 
our lack of understanding of the solutions concentrating in a point xi where {xi,X2,t) is a couple of 
focus at distance t. We know that the solution reconcentrates after a time t in the other focus X2 but 
we do not know precisely how : can it split into several concentrating waves on X2 with different "rate 
of concentration"? That is to say with some different x„ converging to X2 but which are orthogonal. 

Before getting into the proof of the proposition, we state two lemmas that will be useful in the 
proof. Using the notation of Definition II. H denote 

S^{v) = sup I ||V(/9||^2('p^ ,D\vh f, up to a subsequencej 

where the supremum is taken over all the sequences x in M. 
If Vh eL°°{[-T,T],H^{M)), we denote 

5{v) = sup|||V93||^2(7;^^^^) ,Dlvh{th) (f, up to a subsequencej = sup S"" {v{th, ■)) 

where the supremum is taken over all the sequences x = (x/j) in M and t = (t^) in [—T,T]. 
Lemma 2.5. Let ^ G C°^{M). Then, there exists C > such that for any v, we have the estimate 

The proof is left to the reader. 

Lemma 2.6. There exists C > such that for any v = (vh) a bounded strictly (hn) -oscillatory 
sequence in H^{M) 



Urn \\vhhe<C6^v)'/^ Ijm \\vh\\]iLi) ■ 



Proof. This lemma is already known in the case of where the definition of 6^3 is the same except 
that Dj^ is only considered in the trivial coordinate chart. It is estimate (4.19) of |20j in the case of 
a 1-oscillatory sequence, which can be easily extended to (/in)- oscillatory sequence by dilation. 
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Let G C^{Ui) associated to a coordinate patch By Proposition II. 2| ^I'^iVh is still [hn] 
oscillatory and we can apply the estimate on M^. We get 



lim Wi^iVh\ 

n— >+oo 



<C5l^{^,,^ivf/^ lim 11$**, 



n— >-+oo 



Vh 



1 1/6 



< C6^4'^*^,v)'/' lim 



I ll^/*^ 



Then, by definition of the convergence D^, we easily get 
We conclude using Lemma 12.51 and partition of unity. 



□ 



Lemma 2.7. Let T > 0. There exists C > such that for any sequence v = (vh) (hn) -oscillatory, 
solution of the damped linear Klein- Gordon equation on M with bounded energy, we have 



lim \\vh\\LoonrpT] LHM) ^ C^i^)^^'^ '^^^ Wivhi^) , dtVh{0))\\ 



1/6 



Proof. Let th be an arbitrary sequence in [—T,T]. We apply Lemma [2.61 to the sequence Vh{th) and 
get 

li^ \\vh(.tH,-)\\L'^ < C6^v{th,-))'^' li^ \\Vh{th)\\]itM) 



< C6{vf/\\mi\\{vj,i0),dtVh{0m£ 



1/6 



n— >+oo 



by definition of 5 and by energy estimates. 



□ 



Proof of Proposition \2.^ It is based on the same extraction argument as f2] and [17] : the concentra- 
tion will be tracked using our tool and we will extract concentrating waves so that 5{v) decreases. 
We conclude with Lemma [2.71 to estimate the L°°{L^) norm of the remainder term. 

More precisely, if 6{v) = 0, Lemma [2.71 shows that there is nothing to be proved. Otherwise, pick 

{xf^\tf^'^) converging to {x"^ ,t^^) and {(p^'^\'il>^^^) G S^^, such that 



and 



> 



2 1 

> -S(v) 



h-¥0 



The existence of the weak limit ip^^^ (up to a subsequence) is ensured by the boundedness in 
of dtVh (considered in a coordinate chart) by conservation of energy. 

Then, we choose p^^^^ as the damped linear concentrating profile associated with 
[((/j(i),^(i)),/i,x(i),t(i)] (actually, we pick one representant in the equivalence class modulo sequences 
converging to in the energy space as in Definition I0.2p . Remark here that the assumption t^j^^ G 

[— T, T] ensures t^ G [— T, T], which will always be the case for all the concentrating waves we 
consider. Then, we give a lemma that will be the main step to the orthogonality of energies. 



Lemma 2.8. Let w 



(1) 

h 



Vh-Ph ■ Then, 



\\iv,,dtv,)it)r,= (pl'\a,p«)(t) 



+ 



where the o(l) is uniform for t in bounded intervals. 
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Proof. We first compute the energy at time t^j^\ We denote B the bihnear form associated with the 
energy : 



Bia,b)= [ , 
Jm 



ab + Va-Vb + dta dtb 

Jm 

We have to prove 

By weak convergence to in of Vh, p^f^^ and we can omit the term a b of B. By 

construction and Lemma 11.61 we have D9'\vh,dtVh){t9'^) {ip^^^i/j^^^) and D^f^\p^^\ dtp^}^) 

(93(1), ^(1)). Therefore, D^^\w^f^\dtw'f^^){t''f^^) (0,0). Lemma [IH] gives the expected result. 

Remark that if a = 0, this is just a consequence of the conservation of scalar product for solution of 
linear wave equation. □ 

We get the expansion of Uh announced in Proposition 12. 21 bv induction iterating the same process. 
Let us assume that 

Vh{t,x) = Y!j = iPn\t,x) +w'i\t,x), 



(33) \\{vh,dtVh)\\l = Y.\=i iPh\dtp'i^) 1+ {w1\dtw1^) I + o{l), uniformly int, as h^O 



e s 

and where p^^'^ is a linear damped concentrating wave, associated with data [{^'^^\4'^^'')TlLi^^\tL'^^] 
mutually orthogonal. 

We argue as before : we can assume 5{w^^^) > and we can pick ((^('+-^), ^('+^)), x*^'+^), such 
that : 



(34) 



(0 



L2(T (,+i)M) 



+ 



(0 



1 



L2(T(,+i)M) 2 



(4^ 5,4'^ ) (tr (^('+^) , ^('+^) ) . 



and we define pll~^^^ as a linear damped concentrating wave, associated with data 
[((^('+1)^ ^('+1)^^0+1)], Again, Lemma [23] applied to w^^^ and pf^^^ implies estimates ([32 



with w^^^^^ = wf' — pf'*'"'^'*. 



Let us now deal with estimate (131 p . Lemma 11.41 combined with energy estimates gives for some 
C > only depending on T and a 



+ 



From this and estimate (j32p . we infer 



L2(T(,)M) 



< c 



j=l \ 

So, the series of general term ^HVc/^^-') 



(i) 



{pf,dtpf)t=,\^ + o{l). 



<C\ini\\[uh,dtUh)\\l<C. 

h—^O 



\Iht ,,)Af) + Ik^'^Wl^T (,)M) ) converges. Using estimate 



34|) . we get 



lim(5(ii;W) = 0. 
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Lemma 12.71 yields 



lim 



w 



(0 



0. 



L°°([-r,T],L6(M) Z-5>oo 

This completes the proof of the first part of Proposition 12.21 Let us now deal with the orthogonality 
result. We will need the following two lemmas. 

Lemma 2.9. Let {x^^\t^^'^) {•!±'^\t^'^^)- Letvh be an h-oscillatory sequence solution of the damped 
linear wave equation such that 

(35) Z)«K,a,t;,)(^)) - V'(^)). 
Then, there exists {(p^'^\ip^'^^) such that, up to a subsequence 

(36) K,att;,)(4'V^Jv'^'^V'('^)- 



Moreover, we have 
(37) 



Proof. First, we assume x^^^ = x^'^\ By translation in time, we can assume t^^-* = 0. The non 

(2) 

orthogonality assumption yields, up to extraction, t\ /h = C + o(l) with C constant. 

Let il)) £ £oo arbitrary and p^ the linear damped concentrating wave associated with 
[{'f,'4^)ih,^^\^]. We use the equivalent definition stated in Lemma [T3] : ([35]) is equivalent to 



V^;/,(0) • Vm(0) 



M 



M 



dtvhmtPh{^) 



(2) 

As both Vfi and ph are solutions of the damped wave equation on M and t, — > 0, we have by 



equicontinuity (see Lemma ll.9p . 



v^,(tf )-vp,(C)+ / dMt):')dtPH{t):') 



.(2) 



M 



.(2)^ 



f(2)^ 



M 



Let V, w satisfying on T^^M 
Conservation of the scalar product yields 



V'(^)'0= / Vv{C)-Vw{C)+ [ dtv{C)dtw{C). 



But according to CoroUarv 12.11 (ph, dtPh)i,_A2) is a concentrating data according to 

[{vu{C),dtw{C)),h,x^^^]. Since the wave equation is reversible and {(f, ip) is arbitrary, we have proved 
that for any concentrating data {fh,gh) associated with [{'f,'4),h,xj^^], we have 



Vvh{t'i'>)-Vfh+ I dtvh{t'^')g, 



IM 



M 



Vv{C)-V^ + 



dtviC)i^. 
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„(1) 



.(2) 



by taking ((^(2)^^(2)^ ^ {v{C),dtv{C)) which satisfies dST]) by 



This gives the result for 
conservation of the energy. 

In the general case x^^'^ we have in a local coordinate chart and up to a subsequence 



.(2) 



- x^^^ + {D+o{l))h where D G T^^M is a constant vector. We remark that if a bounded sequence 
satisfies D\l\fh,gh) (v',V'), it also fulfills I)f + V(. + ^)). □ 



?t-s>0 



We will also need the following lemma which is the analog of Lemma 3.7 of [17]. We keep the 
notation of the algorithm of extraction for further use. 

Lemma 2.10. Let {j, /} G {1, • • • ^K}^ be such that 



Then, Df^'-\wf^^\dtwl^^''){i}C^'') - ^mp/^es D)l\wX^^'\dtwX^^''){tl') - 0. 



{X + l)w,(E' + l)^ 



i(i)r.„(^+i) 



Moreover, if we assume 



(i) _ Af) 

OO t'OO 



focus 



(see Definition \UJ^, then D^^\p^l'\dtp)[ '){t 



(0,0) for any concentrating wave p)^ ^ associated with [{^p^^ \ ^)^hL^-!±'^ Kt!'^ ■*]• 



Proof. The first result is a particular case of Lemma 12.91 The proof of the second part is very similar 
to Lemma 3.7 of [17]. To simplify the notation, we can assume by translation in time that t^j^ ^ = 0. 
We have to distinguish two cases : time and space orthogonality. 



In the case of time orthogonality, that is 



+00, we first prove D^^^^\p^^ ^)i'^h^) 







(recall that the exponent 1 in -D^''-"'^ means that we only consider the part of the weak limit). 
Thanks to the nonfocusing assumption. Lemma 12.31 yields 



We choose {U,^u) some local chart around x'^ and G Cq°{U) equals to 1 around x'^' . Then, 

(here, we have identi- 



(i) 



and /i2 



^uPh'\t'f!\xh + hx) 

fied '^uPh ^ with its local representation in M^). In particular h^'^up^^ \t^il\xh + hx) and 
^h^''\Ph ^)i'ih^) ~^ 0- Now, we want to prove more precisely dI^\pI^ \dtP^h ^ 0. Suppose 

d\1\p^P ,dtp^P){t^l'') ^ (0,^). Take s G M arbitrary. = + «^ fulfills the same assump- 

+00 and the nonfocusing property \too\ < T focus- So, we conclude similarly that 



tion 



(Ph )i^h ) 0- But the proof of Lemma [2J] gives that D]^>ip)^ >,dtp)^ Wii') ^ {v,dtv){s) 
where v is solution of 



□oo^^ = 0, {v,dtv){<d) = {<d,^) 

So, we have v{s) = for any s G M, which gives -0 = and D^\p^f^ \ dtp^^ ^)(4^'') ^ 

In the case of and space orthogonality, Lemma 12.91 allows us to assume that tf^ 

0- lu local coordinates, we have 



X — x 



if) ' 
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If a^oo^ 7^ , the conclusion is obvious. If it is not the case, take g G Cg°(M^). For the first part, 
we have to estimate 

/ Jj) _ Jf) \ 
+ U+ ^ ^ " \g{y)dy 

which goes to as /i tends to because g is compactly supported. The same result holds for the 
second part for dtp^/^ ^ . □ 

Let us come back to the proof of the orthogonality of cores in Proposition 12.21 Define : 

=max{j G {l,---K}\ A 

assuming that such an index exists. 

We list a few consequences of our algorithm 

(38) (vp('+i),V^'+^))with(^('+i)/Oif/<i^ 

(39) 

(40) 4-) = y: ^'i^'+^r'^ 

The definition oi pf^ and Lemma 11.61 implies D^^^ {pj^ ,dtp'y^){t^l^) {ip^^\ip^^'^). Then, we get from 
(f38|) and 1^ that d['+^^(w['+^\ att(;{^'+^^)(tiJ+^^) ^ (0,0). We apply this to / + 1 = ji^ and it gives 
D[^+^^(u;Sf'^\(9tu;Sf'^^)(tJ/^+^^) ^ (0,0) thanks to the first part of Lemma [2T0] and the definition of 
JK- 

The definition of Jk and the second part of Lemma 12.101 gives Dj^^^~^^\pf^\ dtpl[^){t^j^^^^ ) 
(0, 0) for 3K + l<l<K. 

To conclude, we "apply" to equality (00]) and get D];^^'^^K^^''\tf'^^) 99(^+1) while 

we have just proved D^f^'^^\w^^'^\dtW^^'^^){t^^'^^^) (0,0) which is a contradiction and complete 
the proof of the proposition for 2T < T focus ■ 

In the case of and large times, the orthogonality result is a consequence of the orthogonality 
in short times and the almost periodicity. Denote 

JK = max {i G {1, • • • K}| 3m G Z s.t. (^^^ + mvr, {-iTx^^) U (if 4'""''^)} ■ 

Then, for any + 1 < j < K, we can find m^^^ G Z such that 



< 7r/2 < T focus 
(tf +m(^-)vr,(-l)-<^\(^-)) ±, {tf-''\xf^\ 

and we denote m^^^') G Z such that (tjf'^^ + m(j'«)7r, (-l)™^'^'^!;''^^) (tf xf "*"^^). 
We remark that p^^^ + ui^^'^tt, is still a non zero concentrating data associated with 

[(—1)™*^^ (99, ■(/')((— l)™'''' •),/!, ( — 1)%] thanks to Lemma [2^ (note that it is at this stage that we 
use M = and a = : it is the only case where we are able to describe this phenomenon of 
reconcentration) . So, we are in the same situation as before, and we get a contradition. 

This completes the proof of Proposition 12.21 □ 
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Proof of Theorem \0.3[ We only have to combine the both decompositions we made. Denote Vn 
(and the rest pn^ the hn^ oscihatory component obtained by decomposition (j24p and pn'"^ the 
concentrating waves obtained from decomposition (|30p (and the rest Wn'^^^). We enumerate them 
by the bijection a from into N defined by 



a) < a{k, (3) if j + a < A; + /? or j + a = A; + /9 and j < k. 
For / and Aj fixed, 1 < j < /, the rest can be written 



Let e > 0. To get the result, it suffices to prove that for /g large enough, 



< e for 



all (/, j4i, • • • , Ai) satisfying I > Iq and a{j, Aj) > a{lo, 1) . 

([9]) can easily be deduced from the same orthogonality result in the both other decomposition. 



In particular, it gives that the series of general term lim {pn'°^\ dtPn'"^^ 

In particular, we can find Iq large enough such that we have 



)t=o 



is convergent. 



(41) 



J2 lim ipii'''\dJr^)t=o 



< e. 



o-(j,a)>(T(/o,l) 

Moreover, for Iq large enough, we have for I > Iq 



lim 

n— >oo 



< e. 



L°°(L6) 

Then, for any I >Iq, one can find one Bi such that for any 1 < j < I, Aj > Bi implies 



lim 



< e/l. 



The rest can be decomposed by 



where 



q{j,Al,- Al)) 



l<j<lAj<Bi 



j=l Aj<a<Bi 



Since S'n ' ' is solution of the damped wave equation, energy estimates and Sobolev embedding 
give 



lim 



c{j,Ai,-,Ai)) 



L°°{L^) 



< C lim (5^-''^^'-'^'»,5tSi-''^i' 



^ ^E E 

j=l Aj<a<Bi 



t=0 



2 

£ ' 



t=Q 
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where we have used ahnost orthogonaUty in the last estimate. But the sum is restricted to some 
(j, a) satisfying a{j,a) > a{j,aj) > a{lo,l) and is indeed smaller than Ce thanks to (|^T]) . 



Combining our estimates, we get that lim 



{l,Ai,-,Ai) 



is smaller than (2 + C)e for all 



L°°(LS) 

{l,Ai, ■ ■ ■ ,Ai) satisfying / > Iq and a{j,Aj) > (j{lo, 1). We get the same estimates with the L^(L^^) 
norm by interpolation between L°°{L^) and L'*(L^^). The second norm being bounded by Strichartz 



estimates and the fact that Wn 



is uniformly bounded in the energy space. 



□ 



We also state a few consequences of the algorithm of Theorem 10.31 that will be used below. The 
following both lemmas use the notation and the assumptions of Theorem 10.31 . 

Lemma 2.11. Let 2T < Tjocus- For any I G N and 1 < j < I, we have, with the notation and 
assumptions of Theorem \0.3\ 

z)i^-)rf,5rf)(4^-))-(o,o). 

Proof. Assume D^J^ (w^^ , dtw'n^ ) {tl^^ ] 
to write for L > I 



'i/')- We directly use the decomposition of Theorem 10.3 



w. 



(0 



i=l+l 



In case of scale orthogonality of hn^ and hn\ for /+1 < i < L, we have directly Dn\pn \ dtPn^){tn^) 
(0,0). Otherwise, if /i^f'^ = h'^^ and {x^^\t^^^) ±h {x^'Kt^^), LemmaElO] gives the same result. There- 
fore, Dn\w'^\dtw^n'^){tn^) i^,'^)- Since lim — > 0, we have (/? = 0. We 

n-!>oo L°°{[-T,T],L'^) L->oo 

finish the proof as in Lemma 12.101 We use the same argument for times tn'^ + shn'^ and get if) = Q 
by the proof of Lemma 12.91 Remark that Lemma 12.91 requires that w'^ is strictly /i^f^ -oscillatory, 
but this can be easily avoided by decomposing Wn"* = fn + 9n with /„ (/ii^^)-oscillatory and gn 
(/in^) -singular. 

Lemma 2.12. With the notation and assumptions of Theorem \0.3\. we have, for any j G N 



□ 



lim 

n— >oo 



P'. 



(i) 



- '^JSlll^"lli^^([-r,T],LiO) 



L5([-T,T],LiO) 
where C only depends on the manifold M . 

Proof. We first assume 2T < T focus- Actually, in the case of M^, the result is proved using the fact 

that the ph, are some concentration of some weak limit of a dilation of Vn ■ The proof for a manifold 
follows the same path with a little more care due to the fact that dilation only have a local meaning. 
For any e > 0, we prove 



lim 



< C lim ||^'n||L5cr_T'ri rioi + 



L5([-T,T],LiO) 

We use the decomposition of Theorem 10.31 and choose / > j large enough such that 



lim 

n~^oo 



L5{[~T,T],L^0) 



< £. 



Let ^ij be a cut off function related to local charts (U,^u) such that ^ui^) 
^u{x) = around any ^ xio- 



1 around x-L, and 
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For each 1 < i < Z, we decompose [— T, T] = I^'f U I'^'f U I^'f according to Lemma 
For any i such that = x^, for A large enough, we have 



(42) 



hm 

n— >oo 



Moreover, Lemma [2.21 yields for A large enough 



< e/l. 



(43) 



lim 



< e/l 



where v^^\t,x) = ^i=$^^c;(x)f (*) 



j on a coordinate patch and v^"^' solution of 



(44) 



U^j = on M x T M 

(?;»(0),5tw(*'^(0)) = ((^«,'0(*)). 



Thanks to (1421) and (j43p , the conclusion of the lemma will be obtained if we prove 



,(i) 



L5(1R,L10(T j M)) n^oo 



< lim ll^nll rSi-r.T^Tii riO'i + Ce. 



We argue by duality. Take / € C^(M x Tj M) with 



L5/4(R,LiO/9) 



1. 



From now on, we work in local coordinates around xsb and we will not distinguish a function 
defined onU d M with its representant in w T j M. Denote the operator defined on functions 
on Mj X by 

W^g{s, y) := \[hig{ti + Ks, t{ + Ks). 

(i) 

The definition of vii in local coordinates yields 
On the other hand 



xR3 



{W^^ull-T,T]P'n)f = / 



f. 



For any 1 < i <l, with x^ 7^ xso , using again Lemma [2.31 and 1 2.21 and the fact that we can choose 
with ^[/(^^S) = 0, we easily get 



lim 

n—>-co 



L5([-T,r],LiO) 



0. 



So for n large enough 



< (||wn|lL5([_r,T],Lio) + 2e ) + 



C/^[-T,T] 



E 
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But for i ^ j, = using (j42p and then (j43p . we have for A and n large enough 



< 



< 







/ 












/ 









+ e/l 
+ 2e/L 



These terms are actually 

/mxR3 

^l{hix + xi)l , 



(0 



/r 



Since this expression is uniformly continuous in v'^ E L^(M, L^^(M^)), we may assume v'^ in Cq^(Rx 
Then, if ^ — > 0, we have 







f 


JRxIR^ 


n,i 





If rf- — > oo, the change of variable s 

/ 



V 



y 



f 



0( 



hi 
hi 



f,^ gives 

""n 

-7/2 



)• 



If hn = the space or time orthogonality yields that the integral is zero for n large enough. 
In conclusion, for any / £ C^(IR x M'^) with ||/||i5/4(]g ^io/9) = 1, we have proved : 



This gives the expected result by duality. 

The case of is proved by considering subintervals of length smaller than Tjocus where the 
former result can be applied. □ 



2.2 Nonlinear profile decomposition 

2.2.1 Behavior of nonlinear concentrating waves (after S. Ibrahim) 

In this subsection, we recall the description of nonlinear concentrating waves. As explained in the 
introduction, the behavior for times close to concentration is ruUed by the scattering operator on 
with a flat metric. So, we first state the existence of the wave operator on M^, following the notation 
of [2]. We state it for any constant metric on the tangent plane Tx^M w M^. 

Proposition 2.3 (Scattering operators on M^). Let x^o G M and the d'Alembertian operator 
(constant) on T^^M w induced by the metric on M. To every solution of 

r D^v = on RxTx^M 

\ {v{0),dtv{0)) = {ip,i;)e£x^. 
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there exists a unique strong solution u± of 



j DooUi = -\u±\^u± on M X T^^M 
] lim \\iv-u±,dtiv-u±))it)\\, =0. 

The wave operators 

: {v,dtv)t=o H> iu±,dtu±)t=o 

are bijective from onto itself. 

The scattering operator S is defined as S = o 

The analysis of nonlinear concentrating waves computed by S. Ibrahim in |21j shows that there are 
three different periods to be considered : before, during and after the time of concentration. Roughly 
speaking, for times close to the concentrating time, the solution is closed to nonlinear concentrating 
waves on with flat metric and without damping, as described in Bahouri-Gerard [2]: in the fast 
time hnt, it follows the scattering on R^. Before and after the time of concentration, the nonlinear 
concentrating wave is "close" to some linear damped concentrating waves as defined in Table 12.11 
below. This is precised in the following theorem whose proof can be found in S. Ibrahim [21j. Yet, 
in |21| , the result is stated for an equation without damping and we give a sketch of the proof in the 
damped case in Section [2.2.21 

Theorem 2.1. Let v = [{^,'4^),h,x,t)] be a linear damped concentrating wave. We denote by u 
its associated nonlinear damped concentrating wave (same data at t = 0). There exist three linear 
damped concentrating waves denoted by [{^i,ipi),h,x,t)], i = 1,2 or 3 such that : for all interval 
[—T,T] containing t^o, satisfying the following non-focusing property (see Definition \0.1\) 

(45) mes {F^,x^,s) = Vx G M and s / such that too + s G [-T, T] 
we have 

(46) lim|||nn - [((/7l,V'l),Zil,:^,*)]|||rl>A ^ — > 

(47) lim|||ttn - [((/93,V'3),Zl,^,t)]|||r3,A > 

n J-n A— s>+oo 

where, ll;^ = [-T, tn - Ahn] and In'^ =]t„ + A/i„, T]. 

Moreover, for times close to concentration /„' = [t„ — Ahn,tn + A/i„], we have 

(48) VA > 0, lim|||u„ - u;n||L2,A = 

where Wn{t,x) = ^ij(x)-^=w ^ on a coordinate patch and w solution of 



(49) 



UooW = -\w\'^w on M X T^^M 

{W{^),dtwm = (<P2,V'2). 



where 'd^ corresponds to the frozen metric on T^^M . 

The different functions {(pi,^i) are defined according to Table [2n\ following the notation of Propo- 
sition\KR 
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lim ^ 








— oo 

















oo 









Table 1: Transformation of the profile through a focus 



Remark 2.2. Note that the transition from the first column to the third one represents the modifi- 
cation of profile due to the concentration and the concentrating functions are modified according to 
the scattering operator S. To go from the first column to the second one, we apply the operator 
while we apply ilT^ to get from the second to the third one. 



Remark 2.3. The behavior for times close to concentration is not written this way in the article ]21^ 
of S. Ibrahim, hut is a byproduct of its proof. We refer to the next section which contains a sketch of 
the proof. 

Corollary 2.2. A nonlinear damped concentrating wave is strictly (h)-oscillatory with respect to 
Am and hounded in all Strichartz norms, uniformly on any bounded interval. 

Proof of Corollary \2.2l The boundedness of all Strichartz norms is a counsequence of Duhamel for- 
mula and Strichartz estimates once the result is known in the case of L^L^^. On the intervals In^ 

3 A. 

and /„' when is closed to a linear concentrating wave, the result follows from Proposition 11.41 

2 A 

and linear Strichartz estimates. On /„' , qh behaves like a concentration of a nonlinear solution on 
Tx^M. The strict (/i)-oscillation is obvious and the Strichartz estimates follow from global estimates 
on M^. □ 

In the case of S^, thanks to a better knowledge of the behavior of nonlinear concentrating waves 
we can avoid assumption (I45p . This is Theorem 1.8 from |21) . It will allow us to perform the profile 
decomposition for large times. 

Theorem 2.2. Letv = [{f,ip),h,x,t)] be a linear (not damped, that is a{x) = 0) concentrating wave 
on S^. We denote by u its nonlinear associated concentrating wave (same data att = 0). We assume 
that too G]0, 7r[. Then, for all j € Z, we have 



lim 



ln-[S^'^S{ip,ij),h,i-iyx,t)] 



— > 

]t„+jn+Ah„,t„ + {j+l)-!T-Ah„] A->+oo 



where, S = S o A, S^^^ = SoSo---oS,j times and A{ip, ip){x) = —(93, ^){—x). 

Moreover, the cases too G] — '''") 0[ and too = can be deduced similarly to Theorem 12.11 with some 
changes on the concentration data in the same spirit as Table 12.11 



2.2.2 Modification of the proof of S. Ibrahim for Theorem 12.11 in the case of damped 
equation 

In this subsection, we give some sketch of proof for the behavior of nonlinear damped concentrating 
waves announced in subsection 12.2.11 These results are proved in |21| in the undamped case a(x) = 
and so we only briefly emphasize the main necessary modifications of proof. To simplify, we only 
treat the case |^ — > 00. 

71--> + 00 
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Sketch of the proof of estimate J^gp of Theorem \2.1\ : Behavior before concentration. The proof is ex- 
actly the same as Corollary 3.2 of [21j. Wn = Un — Vn is solution of 

OWn + 'Wn + a{x)dtWn = -\Wn + Vn^^iWn + Vn) On In^ X M 
{Wn,dtWn)\t=0 = (0,0). 

Using Strichartz and energy estimates, we are able to use a bootstrap argument if lim llfnlLs/ri.A . iqn 

n—>-oo ^ U'l ) 

is small enough. This can be achieved thanks to Lemma 12.31 and gives the result. □ 

Sketch of the proof of estimate ( (^<g[ j of Theorem \2.1\ : Behavior for times close to concentration. By 
definition of Vn and finite propagation speed, the main energy part of f „ is concentrated for 
times close to too- By estimate (j46p . it is also the case for Therefore, for times t £\tn— A/in, + 
we can neglect the energy outside of a fixed open set and work in local coordinates. Moreover, 
in that case, we can use the norm HHH/xrs instead of |||-|||j and use the fact that is is invariant by 
translation and scaling up to a modification of the interval of time. 

Denote Un (resp Vn) the rescaled function associated to ti„ (resp Vn)^ so that 
Un(t,x) = —^Un I ^T^, "^"T^" 1 . We need to prove lim lllttn — u'llL a Aiyms — > where w is 

V"n \ rin tin J n^OO I il-,i\-JAJR \^oO 

solution of 

OooW = —\w\^w on M X 

{w,dtw)\t=o = (¥'2,V'2) = f^-lv?,^)- 

By definition of ri_, w satisfies \\{w — v,dt{w — v)){t)\\TTiy.T2 — > where v is solution of 

t— oo 

.^Q. r noo^^ = on MxM3 

I {v,dtv)\t=o = (<^,V')- 

Moreover, it is known that Q-{(p, ip) = lim U{—s)Uo{s){ip, ip) where U and Uq are the nonlinear 

s— oo 

and linear flow map. More precisely, by Lemma 3.4 of [21], we have || jtUA — it'll |f_A aIxR^ — ^ where 

' ' ' A— >co 

w\ is the smooth solution of 

OooWa + \wa\^wa = on [-A,A]xM3 
{wA,dtWA)\t=_A = XA{v,dtv)\t=-A- 



where XA is an appropriate familly of smoothing operator. So, we are left to prove 

lim |||n„-?i;A|||[_A,A]xK3 — ^ 0. 

n— >oo ^ ' J A— >oo 

r,A 



We introduce the auxiliary family of functions solution of 

^nUn + h^Un + \Un\^Un = —hna{hnX + Xn)dtUn OU [— A, A] X 
{u^,dtUn)\t=-A = {Vn,dtVn)\t=-A- 

where we have denoted the dilation of the operator □. So it can be written = df — 
Y^i j g^-' {hnX + Xn)dij + hnV{hnX + Xn) ■ V where y is a smooth vector field (note that it is only 
defined in an open set of size 0{h~^) but it is also the case for Un, '^n, we omit the details). 

The proof is complete if we prove 

and 

(52) hm \\\un - Un\\\,^ M^^3 0. 
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We begin with (|5ip . r„_A = n„ — w is solution of 

+ hna{hnX + Xn)dtrn,A = \wA\^WA-\rn,A + WA\^irn,A + WA) 
< -hlwA - hna{hnX + Xn)dtWA + (Doo - ^n)wA 

{rn,A,dtrnA)\t=-A = {Vn - XAV, dt{Vn - XAv))\t=-A- 

A quick scaling analysis easily yields that the operator □„ + /i^ + hna{hnX + Xn)dt satisfies the same 
Strichartz and energy estimates as □ + 1 + a{x)dt for some times of order A. Moreover, following the 
same argument as Lemma 2.1 of [21], we get that for fixed A 

Jirn^ \\-hlwA - hna{hnX + Xn)dtWA + (Doo - l^n)'^^A|Li([_A,A].L2) = 0- 

Thanks to Lemma [2.2| we know that lim \\(vn — XAV^dAvn — XAv)){—A)\\j^^ ,3 can be made arbi- 
trary small for large A. Strichartz and energy estimates give for any rj > —A 

III'^"-A|II[-A,^]xM3 ^ \\{Vn-XAV,dt{in-XAv)){-A)\\^^^^^ 

+ \\-hl^WA - hna{hnX + Xn)dtWA + (^00 - l^n)^^A||iiQ_A,^],z,2) 

I II ||5 , II II II ||4 

+ irn,A|lL5([_A,,,],LiO) + \VnM L'^([-A,ri\,L^'i) II'"^A|Il5([_a,,,],L10) • 

If ll''^A|li5([_A r?] Lio) small enough, a bootstrap gives (jSip on [— A,r/]. We can iterate the process 
by dividing [—A, A] in a finite number of intervals where the bootstrap can be performed. 

For (I52p . we observe that and Un are solutions of the same equation but with different initial 
data which satisfy thanks to estimate (j46p 

lim \\{u^ - Un,dt{Un - Un)){-A)\\c- = lim \\{Vn-Un,dt{Vn-Un)){-A)\\£ . > 0. 

n— >oo " n— >oo A— >oo 

Then, Strichartz and energy estimates allow us to use a boot strap argument on subintervals I such 
that ll^^n ll/^sfT" j;^io-) is small. (I5ip allows to complete the proof. □ 



2.2.3 Proof of the decomposition 



This subsection is devoted to the proof of Theorem [0? 
Let us define the function (3 in the following way : 

Proposition 2.4. Let < 2T < Tf 

ocus (see Definition [O^jp. Let Pn\ 1 ^ J ^ ^, linear damped 
concentrating waves, associated with data [{'f^-^\iJ^-^'^),h^-'\x''-^\t^-''^] (we can have h^^^ = 1 for one of 
it), which are orthogonal according to Definition and such that t'^ G [—T,T]. Denote Qn^ the 
associated nonlinear damped concentrating waves (same data at t = 0). 
Then, we have 



(53) 



lim 

n— >oo 




Li([-T,T],L2) 



Proof. We follow closely Lemma 4.2 of [17J. 




Li([-T,T],L2) l<Ji.-J5</ 



/c=l 



n 



L^([-T,T],L'^) 
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(ik) 

where at least two qn are different. In the case of othogonahty of scales, we use Holder inequality 



k=l 



,(ifc) 



< C \\QnQn\\L°^{[_T,T],L3) 11 ||^' 
Li([-T,r],L2) fc=3 



Uk) 



L3{[-T,T],L^S) 



Then, Corollary 12.21 and Lemma 11.81 yield the result (note that L^L^^ is a pair of Strichartz norm). 
So now, we can assume h\ = h'^ = hn- By Holder and Corollary 12. 2^ we get 



n 

k=l 



rOfc) 



Li([-T,T],L2) 



A;=3 



L5{[-T,T],LiO) - ^ lkn^nlL5/2{[_T,T],L5) ' 



We apply Theorem 12.11 to q\. We obtain three couples ^'), i = 1,2,3 and split the interval 



—iln ■ We first deal with the interval In'^. Denote — [(i/?!, -^i^)] so that 

||9n9n|L5/2(ji,A^^5) < lkn|L5(7l.A^^iO) < C' 1 1 9n " ^l,n 1 1 i5(7l>A_^iO) + || ^^l,n || rs r r^A 



1,A 



L5(/^''\L10) 



So, combining Theorem 12.11 and Lemma 12.31 yields 



Il5/2(/^-'^,L5) 



A— s>oo 



0. 



The same result holds for In^ and we are led with the interval In^ . In the case of time orthogonality, 
say — > +00, the two intervals [tl — Khn,tl + A/i„] and [ti — khn^ti + A/i„] have empty 

intersection for fixed A and n large enough, which yields the result by the same estimates applied to 
q^, once A is chosen large enough. 

We can now assume, up to a translation in time, that = t^. On Jn^, Theorem 12.11 allows us 

to replace by w\{t,x) = ^\j{x)'w^ ("TT^; "TT"^) ^ coordinate patch where is solution of a 
nonlinear wave equation on the tangent plane T^i^M and similarly for q'^. In the first case of space 
orthogonality, that is ^ x^, the result is obvious on the interval In'^ by taking "^jj and with 
empty intersection. In the case = x^, we are left with the estimate of 



/2 



\wl{t,x)wl{t,x)\ 



1/2 



ds < 



-A,A] 



w\t,x)w\t,x + ^^:^-^) 



5\ 1/2 



ds 



This yields the result in the last case of space orthogonality by approximating and by compactly 
supported functions. □ 

In the case of the sphere, we are able to state the same result without any restriction on the time. 

Corollary 2.3. Let M = and T > (eventually large). We make the same assumptions as 
Proposition \2.4\ except for the time T, with the additional hypothesis : 

[hi^\ ( — l)'"x^*\ t*^*^ + rail] is orthogonal to [h^^\x^^\t^^^] for any m £ Z, and i ^ j. Moreover, we 

assume a(x) = (undamped equation). 

Then, the same conclusion as Proposition \2.4\ is true. 



Proof. We build a covering of the interval [—T, T] with a finite number of intervals of length stricly 
less than Tfocus = vr so that on each of this interval / = [a,/3] and for any 1 < i < I, there exists at 
most one m^*^ G Z such that t'^ + m*^*V € /. Moreover, one can also impose a 
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Therefore, a e]4'^ + (m* - 1)tt + Kh^n\t'i^ + m^^\ — Ahn^] for large fixed A and n large enough. 



Theorem yields 



for a linear concentrating wave Vn^ 



£ n— >oo 

[5™ Siif'-^ , V^'^ ) , /i^^^ , x(^) , t )] . In each interval, we are in the same situation as in Proposition [2 
which yields the desired result. □ 

Now, we are ready for the proof of the nonlinear profile decomposition. We give it in a quite 
sketchy way since it is very similar to the one of [2] or |17| . First, we obtain it in the particular case 
where the linear solution is small in Strichartz norm. 

Lemma 2.13. There exists 6i > such that if 

lim^WVnW m[-T,T],LiO) < h 

then the conclusion of Theorem \0.4\ is true. 

Proof. The proof is essentially the same as Lemma 4.3 of [2]. We have to estimate the rest rn^ 
solution of 

f Dri') + rH^ + a{x)dtril^ = f]{u) + Y!j=i Piqi^^ ) - P (u + Y!j=i Qn^ + ^i'^ + ^n^) 
I ir^\dtrli\=o = (0,0). 

We conclude as in [2j using Proposition 12.41 and Lemma |2. 121 which is not immediate on a manifold. 
In the case of and a = for large T, we use Corollary 12.31 instead of Proposition 12.41 . □ 

Once the result is obtained when Strichartz norms are small, we divide [—T, T] in a finite number 
of intervals where the Strichartz norms are small enough. This is done in the following lemma. 

Lemma 2.14. Let 2T < Tj 

ocus- Let (5 > and he a sequence in T, T], ij^'^(Af )), such that 

Jim^||g„||^5([_r,T],Lio) < ^- 

(i) 

Fix also I G N and I sequences of nonlinear concentrating wave qfl , j = 1, /. 

Then, for any 6' > 6, there exists L G N such that for any n G N, we have the decomposition of 
[—T, T] in closed intervals In^ 



[-r,r] = U/0-), 



i=l 



such that the sequence 



satisfies on each interval In^ 
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Proof. We first treat the case / = 1. We divide [—T, T] = In^ U In^ U In^ according to Theorem 12.1 
(one of these intervals being possibly empty). Then, a combination of estimate (146 P of Theorem 12.1 
(comparison with linear concentrating wave) and Lemma 12.31 (non reconcentration) gives for A large 
enough 



lim 

n— >oo 



■in 



L-^{iy,LW) 



<6' -6. 



The same result holds for In^ and we are left with the interval In^- Once A is fixed, we can divide 



[—A, A] in a finite number of intervals such that 11^^11^5(7(0. a lio < S — 6' where w is the function 

defined by equation (I49p of Theorem 12. 11 Then, we replace each by In ^'^ obtained by translation 
dilation. We conclude by the approximation (|48p of qil'^ by translation dilation of w on the interval 



2,A 



□ 



Note that the previous lemma also applies for large times on 5^ with a = by doing a first 
decomposition of [— T, T] in a finite number of intervals of length strictly less than vr. 



(0 



< 5i and use 



End of the proof of Theorem \0-4\ in the general case. We choose / G N such that w 

Lemma [2. 141 in order to be able to apply Lemma [2. 131 on each interval In^ . See [2] or, in the different 
context of Schrodinger equation, [25| . □ 



2.3 Applications 

2.3.1 Strichartz estimates and Lipschitz bounds for the nonlinear evolution group 

Proposition 2.5. Let T > be fixed. There exist a non- decreasing function, A : [0,oo[— t- [0,oo[, 
such that any solution of 



(54) 
fulfills 



□n + u + a{x)dtu 

{u{^),dtum 



I l4 

\u\ U 



on 

{uQ,ui) G £ 



-T, T]x M 



Il8([-T,T],L8(m) + ll^llL5([-T,T],LiO(Af) + II^IIl4([_t,T],L12(M) ^ ^(11 (^0, ^^l) lU) ■ 



Proof. The proof is exactly the same as Corollary 2 of [2j. Using Strichartz estimates, it is enough 
to get the result for L^L^^. We argue by contradiction and suppose that there exists a sequence Un 
of strong solutions of equation (I54p satisfying 



sup \\{uo,n,Ul^n)\\£ < +00, 



''n\\L^{[-T,T],L^O{M) 



+00. 



We apply the profile decomposition of Theorem 10.41 to our sequence. We get a contradiction by 
the fact that the L^{[—T,T],L^^{M)) norm of a nonlinear concentrating wave is uniformly bounded 
thanks to Corollary 12.21 This argument works for times 2T < Tf^cus and can be reiterated since 
the nonlinear energy at times T can be bounded with respect to the one at time thanks to almost 
conservation (we can also use energy estimates once we know u is uniformly bounded in L^L^^). □ 

Lemma 2.15. Let Rq > and T > 0. Then, there exists C > such any solution u satisfying 



(55) 



Du + u + a{x)dtu + \u\'^u 
{uiO),dtu{0)) 
\\iuo,ui)\\£ 



< 



on [-T, T] X M 
{uo,ui) e £ 

Rq. 
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fulfills 

\\{u{t),dtum\L^xH-^<C\\iu{0),dtuiOm^,,^^^ ytG[-T,T]. 

Proof. Proposition 12.51 yields a uniform bound for u in L'^([—T,T],L^'^(M)) and so for V = \u\^ in 
L^{[0, T], L^{M)). We prove uniform estimates for some solutions of the linear equation 



(56) 



Uu + u + a{x)dtu = Vu on [-T, T] x M 
{u{ld),dtum = (no,ni) G L2 X H'^ 



where V satisfies ^3(4/)) — ^(-^o)^- The product of functions in L°° (\—T .,T\, L?') and 

Li([-T,r],L3) is in Li([-T, T], L^/S) and so in /.^([-T, T], i^-^) by Sobolev embedding. Standard 
estimates yields 

ll(^'^t^)llL-([0,t],L2xH-i) < <^ 11(^(0), 9i^^(0))||i2x//-i +<^(t+ ||1^|lLi([0,t],L3)) ll(^'^i^)llL-([0,t],L2xH-i) ■ 
We can divide the interval [— T, T] into a finite number of intervals [oj, aj+i]i=i...7v such that C(t + 
ll^llLi(KA+i],i=^{A/))) < V2- N depends only on Rq and T (not on V). 
Then, on each of these intervals, we have 

<2C\\{u{ai),dtu{ 

We obtain the expected result by iteration. The final constant C only depends on Rq and T since it 
is also the case for N . □ 

Corollary 2.4. Let Rq > and T > 0. For any e > 0, there exists 6 > such that any solution u 
satisfying ( 1551) and \\{uo,ui)\\j^2y^jj~i < S satisfies 

||(u(T),5<n(r))||^,,^_, <e. 

We will also need the following lemma which states the local uniform continuity of the flow map. 
Note that it can be proved to be locally Lipschitz with a slightly more complicated argument (see 
Corollary 2 of [17]). We will note need this for our purpose. 

Lemma 2.16. Let Un, Un be two sequences of solutions of 

□n„ + Un + \Un\^Un = Qn Oli [-T,T] X M 

iun,dtUn)t=o = (M„,o,Wn,i) bounded in 
with \\{un,o - Unfl,Un,i " Un,i)\\£ + \\9n " 9n\\ (^[_T,T],L^) ^- Then, wc havc 

Proof, rn = Un — Un is solution of 

Orn + rn + \Un\^Un - \Un\'^Un = Qn - Qn On [-T, T] X M 



{rn,dtrn)t=0 = {Un,0 - Un,0,Un,l - Un,l). 

Using energy and Strichartz estimates, we get 

1 1 l^'nl 1 1 [-T,T] — ^ II (^n,0 ~ Unfl, Un,l — Un,l)\\^ + C ll^n ~ 5n|| Li ([-T,T],L2) 

+C! \\'^n\\L5([-T,T],L^0) (\\Un\\L5([-T,T],L^0) + ll""™ IIl5([-T,T],L10) 

Using Proposition 12.51 we can divide the interval [— T, T] in a finite number of intervals n — 
[ai,„, ttj+i^n], 1 < i < iV, such that C' (^||n„||^5(7. ^^^lo) + ||u„||'^5(7. ^^^lo)^ < 1/2 so that the third 
term can be absorbed. We iterate this estimate times, which gives the result. □ 
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2.3.2 Profile decomposition of the limit energy 

For u solution of the nonhnear wave equation, we denote its nonhnear energy density 
e{u){t,x) = ^ [\dtu{t,x)\'^ + \Vu{t,x)\'^ + |u(t,x)|2] + ^\u{t,x)f. 

For a sequence u„ of solution with initial data bounded in £, the corresponding non linear energy 
density is bounded in L°°{[—T,T],L^) and so in the space of bounded measures on [— T, T] x M. 
This allows to consider, up to a subsequence, its weak* limit. 

The following theorem is the equivalent of Theorem 7 in [8j. It proves that the energy limit 
follows the same profile decomposition as Un- It will be the crucial argument that will allow to use 
microlocal defect measure on each profile and then to apply the linearization argument. 

Theorem 2.3. Assume 2T < T focus ■ 
Let Un be a sequence of solutions of 

UUn + Un + \Un\'^Un = 

with (un,dtUn){0) weakly convergent to in E. 

The nonlinear energy density limit of Un (up to subsequence) reads 

+00 

e{t, x) = ^ e^-'') (t , x) + e/ (t, x) 
i=i 

where e^^^ is the limit energy limit density of Qn ^ (following the notation of Theorem \0.4\ ) and 

er = lim lim e(tt;^'^) 

where both limit are considered up to a subsequence and in the weak * sense. 
In particular, ej can be written 



ef{t,x) = / n{t,x,d^) 



with 

fi{t,x,0 = l^-iGt{x,0) + l^+{G-t{x,0) 

where Gt is the flow map of the vector field -f^|^|^ on S*M, that is the Hamiltonian of the Riemannian 
metric. 

Moreover, e is also the limit of the linear energy density 

eiin{Un){t,x) = ^ [\dtUn{t,x)\'^ + \\7Un{t,x)\'^] . 

Proof. Proposition 12.51 yields ||'Un||i8([_T xM) — G. Then, compact embedding and Lemma 12.151 
yields ||^in|l2,2([_r2^]xAf) — ^ ^ ®° II'"" \l^{[~TT]xM) — ^ ^ interpolation. Therefore, e is the 
limit of b(un,Un), with 

b{f, 9) = dtf{t, x)Wg{t, x) + V/(t, x) ■ W7){t, x) 
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Now, we have to compute the hmit of b{un, Un) using decomposition (jlOp of Theorem 10.41 We set for 
any / G N 



,(0 



and so 



6K,n„) = + + 26(4'), + 26K,r«) - fe(r«,rW). 

converges to zero as / tends to infinity. 



Because of (jlip . hm 

=(0 



So, if we define Cr' = w* hm ( 26(u„,ri')) — 5(r 



Li([-T,T]xM) 



we have 



0. 



TV i-s>oo 

Let (p S Cg°(] — T, T[xM). For fixed I, it remains to estimate 

I 



-T,T[xM 



E 



-T,T[xM 



Since b{qi^^ , Wn'') is bounded in — T, Tf, L-'^), we can assume, up to an error arbitrary smah, that 
is supported in |t < t^^j or |t > t^^j (replace by (1 — 'I')(t)(^ with ^'(t^^) = 1 and H^HiiQ-jTi) 

smah). On each interval, Theorem 12.11 allows to replace qn^ by a linear concentrating wave. Then, 
we combine Lemma [2.111 and Lemma [1.1 II to get the weak convergence to zero of b{sn\wn^) for fixed 

Lemma [2.101 and the orthogonality of the cores of concentration give D^^^p'^ \dtp)ji ^)(i/f^) ^ 

(0,0) for j 7^ j' and p^^ ^ a concentrating wave at rate [h^^'\t^^'\x^^'^]. Then, the same argument as 
before yields 



n- 

So we have proved that for any / E N 

I 

h{un,Un) e = y^e^^^ 



,0') 



+ elf? + €: 



(0 



where e'w is the weak* limit of b{w^ ,w^) and er satisfies 



,(0 



0. ei) is the weak* limit of 



TV /->oo 

a sequence of solutions of the linear wave equation weakly convergent to in energy space. Therefore, 
it has the announced form using the link with microlocal defect measure (see Lemma [1.9p . 

We get the final result by letting I tend to infinity. □ 

Remark 2.4. The fact that |n„|^ is weakly convergent to is false if we consider the limit in D'(M) 
time by time. For example, for a nonlinear concentrating wave with tn = 0, the weak limit in 
T>'(\ — T, T[xM) of \un\^ is of course still zero but the weak limit of |u„|^(t) in D'(M) is zero ift^O 
and a multiple of a Dirac function ift = 0. So the limit in D'(M) of en\t=o ^■^ '^'^^ same as the 
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one ofb{un, Un)\t=o- This comes from the fact that the limit ofb{un, Un){t) is not equicontinuous as a 
function of t while it is the case for the nonlinear energy. Yet, in the proof, we are only interested in 
its limit in the space-time distributional sense which will be continuous. Actually, the discontinuity 
at t = of the limit of b{un,Un){t) can be described explicitely from the scattering operator. At the 
contrary, the fact that the nonlinear energy density e{t) is continuous in time can, in this case, be 
seen as a consequence of the conservation of the nonlinear energy of the scattering operator. 

3 Control and stabilization 

3.1 Weak observability estimates, stabilization 
3.1.1 Why Klein-Gordon and not the wave? 

In this subsection, we prove that the expected observability estimate 

E{u){0) <C [[ \adtu\'^ dtdx. 

JJ[0,T]xM 

does not hold for the nonlinear damped wave equation Du + dtu + u^ = (in the simpler case a = 1), 
even for small data. It explains why we have chosen the Klein-Gordon equation instead. The main 
point is that for small data, the nonlinear solution is close to the linear one which has the constants 
(in space-time) as undamped solutions (which is obviously false for \3u + u = 0). 

We take a = 1 and initial data constant equal to (£,0). The nonlinear wave equation takes the 
form of the following ODE 

f u + ii + u^ = on [0,r] 
\ (n(0),ti(0)) = (£,0). 

Decreasing of energy yields for any t>0 

E{t) = -u^ + -v^{t) < E{0) = -e^ 
\ J 2 6 ^ ^ ~ ^ ^ 6 

and so 

\u{t)\ <£ \/t>0. 

Then, c = u is solution of 

( c + c+u^ = on [0,r] 
\ c(0) = 

Therefore, 

c{t) = - f^e-^'-'^u^s) ds 
and \u{t)\ = \c{t)\ < . 

For any T > 0, we have 




Therefore, the observability estimate 

Te^^ > r \uis)\'^ > CE{0) = cje' 
Jo o 

can not hold if e is taken small enough. 
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3.1.2 Weak observability estimate 

As explained in tlie introduction, the proof of stabilization consists in the analysis of possible se- 
quences contradicting an observability estimate. The first step is to prove that such sequence is 
linearizable in the sense that its behavior is close to solutions of the linear equation. 

Proposition 3.1. Let oj satisfying Assumption\KM o,nd a G C°^{M) satisfying a{x) > rj > for all 

X G w. Let T > Tq and Un be a sequence of solutions of 

( \JUn + Un + \Un\^Un + a{x)'^dtUn = OU [0,T] X M 
1 {Un,dtUn)t=0 = {uo^n,Ul^n) ^ £■ 

satisfying 

iuo,n,ui,n) weakly in £ 



n—^oo 

|2 



(58) / / \a{x)dtUn\ dtdx — > 

JJ[0,T]xM 

Then, Un is linearizable on [0,t] for any t < T — Tq, that is 



where Vn is the solution of 



= on [0, T]x M 

{Vn,dtVn)t=0 = {uO,n,Ul^n)- 



Proof. Denote = sup |s G [0, T] lim — v„|||[q = o| and we have to prove > T — Tq. If 

it is not the case, we can find an interval [t* — e,t^: — e + L] C [0, T] with Tq < L < Tfocus and 
< 2e < L — To (if t* = 0, take the interval [0,L] C [0,T]). Then, Lemma l3.ll below gives that 
Un is linearizable on [t^, — e, + e]. We postpone the proof of Lemma l3.ll and finish the proof of 
the proposition. The definition of t=K gives lim |||nn — Vn|||[ot,-e] = we have proved that 



n— >oo 

lim lllun — I'nlllrt f =0 where Vn is solution of 



UVn = ; {Vn, dtVn)t=U-e = {Un^dfUn) 



t=U 



Since the norm |||-||| controls the energy norm, this easily yields lim |||n„ — ^^n|||[ot,+e] = which is 
a contradiction to the definition of . □ 

Lemma 3.1. With the assumptions of Proposition [Ql Consider the profile decomposition according 
to Theorem \0.4\ of Un on a subinterval [to, tQ + L] C [0, T] with Tq < L < Tfocus- 

Then, for any < e < L — Tq, this decomposition does not contain any non linear concentrating 
wave with t'^ G [to, to + Un is linearizable on [to, to + e]- 

Proof. To simplify the notation, we work on the interval [0, L]. Moreover, since a{x)dtUn tends to 
in L^L^, Lemma [2.161 allows to assume with the same assumptions that n„ is solution of the nonlinear 
equation without damping. Proposition 12.51 and Lemma 12.151 (with Rellich Theorem) give that ti„ is 
bounded in L^([0, T] x M) and convergent to in L'^{[0, T]xM). Therefore, n„ tends to in L^([0, T] x 
M) and so |M„|^ti„ is convergent to in L'^/5([0,r] x M) ^ L'^I'^{[Q,T] x M) ^ H-l{]0,l[xM). 
Then, if we consider the (space-time) microlocal defect measure of n„, the elliptic regularity and 
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the equation verified by Un gives that /x is supported in = ^ in the hnearizable case. So, 

combining this with ([58|) . we get 



OmHU]0,L[ 



xco) 



Using the notation of Theorem 12. 3| this gives e = on ]0,L[xa;. Since ah the measures in the 
decomposition of e are positive, we get the same result for any nonhnear concentrating wave in the 
decomposition of Un, that is 

qij) ^ hi HU]0,L[xu;) 
and if /x(^) is the microlocal defect measure of qn \ we have 



(59) 



/i(j') = in S*{]0,L[xuj). 



Assume that tl^ G [0, e] for one j G N, so that the interval [t^,L] has length greater that 
Tq. Denote the linear concentrating wave approaching qn^ in the interval In^ according to the 
notation of Theorem 12. H so that for any < t < L (here, we use the fact that L < Tjocus)^ we 
have 



In in 



0. 



[t,L] rt-s>oo 

In particular, /u(j') is also attached to pi^ on the time interval ]t^^,L\. Since is solution of 
the linear wave equation, its measure propagates along the hamiltonian flow. Assumption 10.21 and 
L — t'^ > Tq ensure that the geometric control condition is still verified on the interval [t^^ , L] 

which gives fj,^^'^ = when combined with (I59p . That means pn^ = and so Qn^ = as expected. 
Then, for the profile decomposition of n„ on the interval [0, L] (here the weak limit u is necessarily 



zero 



we have proved that tn^ G 



£,L]. Then Theorem 12.11 and L < T focus provides a linear concentrating 



(j) _ U) 

Hn Pn 



wave Pn such that lim 

n~>o< 

Moreover, the conclusion of Theorem 12.11 give lim 

n— >o<; 

lim ||itn|lL5([o e] Lio) = ^ therefore 



fo 

while Lemma 12.31 give lim 



(i) 



0. 



(0 , (0 



L5([o,£],LiO) 

0. This finally yields 



\Un\ ^n||Li([o,e],L2) 



This gives exactly that n„ is hnearizable on [0,e]. 



□ 



We are now ready for the proof of some weak observability estimates. We recall the notation 
E{u) for the nonlinear energy defined in ([2]). 
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Theorem 3.1. Let oj satisfying Assumption 10.^ with Tq and a G C°°{M) satisfying a{x) > r/ > 
for all X G Let T > 2Tq and Rq > 0. Then, there exists C > such that for any u solution of 



' nu + u + \u\'^u + a'^{x)dtu = on [0,T] x M 
{u,dtu)t=o = {uo,ui)££ 
||(no,ni)||^ < Ro 



(60) 
satisfies 

E{u){0)<c{ [[ \aix)dtu\^ dtdx + \\{uo,ui)\\^2^H~^ E{u){0)] . 

\JJ[0,T]xM j 

Proof. We argue by contradiction : we suppose that there exists a sequence ti„ of solutions of (j60p 
such that 

\ 1 

a{x)dtUn\ dtdx +\\{uo^n,Ul^n)\\L2^H-l E{un){0) ] <-S(n„)(0). 
[0,T]xM J ™ 

Denote a„ = (£^(ti„)(0))^/^. By Sobolev embedding for the L^ norm, we have a„ < C{Ro). So, up 
to extraction, we can assume that an — > a > 0. 
We will distinguish two cases : q > and a = 0. 

• First case : an — > a > 

The second part of the estimate gives \\(uq n^ui n)\\ r^s. m-i — ^ ^^'^ (uq n,ui n) in 

X L^. Therefore, we are in position to apply Proposition 13.11 and get that Un is linearizable on 
an interval [0,L] with L > Tq. We get a contradiction to a > by applying the following classical 
linear proposition, which can be easily proved using microlocal defect measure as in Lemma l3.ll 

Proposition 3.2. Let uj satisfying Assumption \0.S\ with Tq. Let T > Tq and Vn be a sequence of 
solutions of 

j Dvn = on [0, T]x M 

1 {vn{0),dtvnm ^ Om£ 

satisfying 

\a{x)dtVn\'^ dtdx — > 

[0,T]xM 

Then, {vn{0) , dtVn{0)) — > for the strong topology of x L^. The same result holds with Dun 
replaced by + Un ■ 

• Second case : a„ — > 

Let us make the change of unknown Wn = Un/o-n- Wn is solution of the system 

(61) + a^{x)dtWn + Wn + an\Wn\'^Wn = 

and 

2 1 

\a{x)dtWn\ dtdx < — . 

[0,T]xM 
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We have for a large constant C > depending on Rq and for all t G [0,T], 



^ \\{Un,dtUn)\\l < E{Un) < C || 5tn„) ||| . 



Therefore, we have 



(62) \\[Wn[G),OtWn[^))\\s = ^ > ^= > 0. 



||(n „(0),a,n„(0))||g ^ 1 
V^(n„(0)) - 

Thus, we have ||(tt;„(0), 9fW„(0))||£- 1 and {wn-,dtWn) is bounded in L°°([0, T], <f ). 
Applying Strichartz estimates to equation (j6ip . we get for C = C{Rq) > 

ll^n|lL5([o,T],LiO) - + ll^nllL5{[0,T],LiO)) 

Then, using a bootstrap argument, we deduce that ||wn||i5([o t] Li") bounded and therefore 

nwn + Wn in Li([0,r],L2). 
Proposition 13.21 yields that Wn converges strongly to some w solution of 

Dw + w = 0; dtw = on w. 

We deduce as in J. Ranch and M. Taylor [32j or C. Bardos, G. Lebeau, J. Ranch [3j that the set of 
such solutions is finite dimensional and admits an eigenvector w for A. By unique continuation for 
second order elliptic operator, we get dtW = 0. Multiplying the equation by w and integrating, we 
obtain w = (note that, at this stage, the choice of the Klein-Gordon equation instead of the wave 
equation is crucial to avoid the constant solutions). We conclude that (zi;„(0), 9tt(;„(0)) tends to 
strongly in £ which gives a contradiction to (I62p . □ 

3.2 Controllability 
3.2.1 Linear control 

In this section, we recall some well known results about linear control theory and HUM method. Let 
(<I>o,<I>i) G X H^^. We solve the system 

+ $ = on [0, T] X M 
{^,dt^)\t=o = (^o,^i)- 



(63) 
and 
(64) 

The HUM operator S from x H^^ to x is defined by 

5($o,^i) = {-dtv{<d),vm. 
Lemma 3.2. If uj satisfies the geometric control Assumption \0.1\ then S is an isomorphism. 



\Jv + v = a^^ on [0, T] x M 
{v,dtv)\t=T = (0,0). 
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Proof. Multiplying equation (|64p by "3>, integrating over [0,r] x M and integrating by part, we get 
the formula 



T 

2 



a$|^ = -/ dtv{0)m + r;(O)9i0(O) = (5($o,^'i),($o,^i)) 
Jm Jm Jm 

where (., .) denotes the duality between x and x H~^. We get the conclusion thanks to the 
following observability estimate which can be proved by the same techniques used in the nonlinear 
problem 



Jo JM 



□ 



3.2.2 Controllability for small data 



Theorem 3.2. Let oj satisfying Assumption \0. 1\ and T > Tq. Then, there exists d > such that for 
any {uq,Ui) and {uo,ui) in x L'^ , with 

\\{uo,ui)\\^ <S; \\iuo,ui)\\£ <5 

there exists g G L°°([0, 2T], L^) supported in [0, 2T] x uj such that the unique strong solution of 

Du + u + lul'^u = g on [0,2T]xM 
{u{0),dtu{0)) = {uo,ui). 

satisfies {u{2T) , dtu{2T)) = {uq,ui) 

Proof. The proof is very similar to except that the critical exponent do not allow to use com- 
pactness argument and we use the classical Picard fixed point instead of Schauder, as done in |9j or 
[28], [29] for NLS. By a compactness argument, we can select a € Cq°{uj) with a{x) > rj > for x in 
u) where Co satisfies Assumption 10.11 Since the equation is reversible, we can assume {uq,ui) = (0,0) 
and take the time T instead of 2T. We seek g of the form a^{x)^ where $ is solution of the free 
wave equation as in linear control theory with initial datum ($o,*^*i) G x H^^. The purpose will 
be to choose the right ($05*^*i) & x to get the expected data. We consider the solutions of 
the two systems 

+ $ = on [0, T] X M 
{^,dt^)\t=o = (^o,^i) 



and 
(65) 



nu + u+\u\'^u = a^^ on [0,T]xM 
{u,dtu\^T = (0,0). 



Let us define the operator 

L : L2 X ^ X L'^ 

^ ($o,^i) ^ L{<^o,^i) = {u,dtu)\t=o- 



We split u = V + ^ with ^ solution of 
(67) 



D^ + if = a2^> on [0,r]xM 
{^,dt^)\t=T = (0,0). 
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This corresponds to the Hnear control, and {—dt^, ^)|j=o = S{^o, $i). As for function v, it is solution 
of 

\Jv + v = -\u\^u on [0,r]xM 
{v,dtv)\t=T = (0,0). 

$ belongs to C([0,r],L2). So, u, v and ^' belong to C([0,r],i7i) n C\[0,T], L^) n L^{[0,T], L^^). 
We can write 

L(^>o, <^i) = K{^o, ^i) + S{%, $i) 

where i^($o, $i) = (-^i^-, ^)|t=o- -^^(^0, ^i) = (-""i, ^^o) is equivalent to ($o, $i) = -5'-^Jr(^>o, $i)+ 
S'~^(— ui, no). Defining the operator B : L'^ x — ^ x ff"^ by 

B{^o, $i) = -S-iK(«>o, $i) + S-H-ui,uo), 

the problem L($05*^*i) = {~ui,uq) is equivalent to finding a fixed point of B. We will prove that if 
||(no,ni)||£; is small enough, i? is a contraction and reproduces a small ball Bji of x H~^. 
Since S" is an isomorphism, we have 

\\B{<^o,'^i)\\ < C{\\K{^o,^i)\\^2y,H^ + \\{uo,Ui)\\^) 

So we are led to estimate ||ii'($0; ^i)llL2x_f/i ~ ||(^' ^t^)|t=o||£-- Energy estimates applied to equation 
(j68p and Holder inequality give 



||(?^,9t^^)|t=o||£- < ^ I|I'"I'*^IIl1{[0,T],L2) - ^ II^IIl5([o,T],L10) • 
But Strichartz estimates applied to equation (|65p give 

II'"IIl5([o,T],L10) ^ C (||a^^|Li([o,T],L2) + ll'"lli5{[0,T],LiO)) ^ ^ (ll (^0, ^l) |Il2xH-i + ll^lli5([0,T],LiO)) 
Using a bootstrap argument, we get that for ||(*J>0; *^'i)|lL2x_f/-i ^ R small enough, we have 

We finally obtain 

ll^('&o,^i)ll <c(\\{<^o,'^i)\\hxH-i + ||(uo,ni)||^) . 



Choosing R small enough and ||(iiO;^i)|lHixL2 — ^/2C, we obtain ||i?(<I>o, ^'i)||£^2x/i--i < R and S 
reproduces the ball Br. Let us now prove that B is contracting. We examine the system 



(70) 



(71) 

We obtain similarly 



□ (u - n) + (n - tt) + |u|^n - |n|^tt = a^{^-^) on [0,T]xM 

{u-u,dt{u-u))\t^'r = (0,0). 

□ (u-u) + (f-u) + |u|^n-|n|''{t = on [0,r]xM 

{v-v,dt{v-v))\t=T = (0,0). 



5($o,^i)-i?(^o,«'i) _ < C||(7;-i),at(t;-5))|t=o| 



< C \\\u\^u — lul^u 



LH[0,T],L2) 

4 , ii~ii4 



(72) < CR^\\u-u 



- ||n - ■u||^5([o,T],LiO) (II^IIl5([0,T],L10) + II'"IIl5{[0,T],L10)) 
L5([0,T],LiO) 
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where we have used estimate ([69|) for the last inequahty. Applying Strichartz estimates to equation 
(1701) . we get 



\u — u\ 



L5([0,T],LiO) 



< C(|||ti| u — Inl^-ulLi 



LH[0,T],L^) 



Li{[0,T],L2) 
($0,^'l) - (^0,^1 



If i? is taken small enough, it yields 



(73) 



u — u 



L5([0,T],L10) 



< c 



($0,^l) - {^0,^1, 



Combining (j72p and (j73p . we finally obtain for i? small enough 



S(«>o,$i)-S(l>o,l'i) 



L'^xH-^ 



L^xH-i 



and S is a contraction for R small enough, which completes the proof of Theorem 
3.2.3 Controllability of high frequency data 



□ 



This subsection is devoted to the proof of the both main theorem of the article : Theorem 10.21 and 

[an 



Proof of Theorem W.SX First, by decreasing of the energy and Sobolev embedding, there exists some 
constant C{Rq) such that the assumption || (uq, wi)||£; < Rq implies 

(74) E{u){t) < C{Ro) and \\{u, dtu){t)\\^ < C{Ro); Vt > 0. 

Fix T such that Theorem 13.11 applies . Then, there exists e > such that for any {uq,ui) satisfying 

(75) \\{uo,ui)\\e < CiRo); ||(no,txi)||^.^^_i < e, 
we have the strong observability estimate 

E{u){0) <C [[ \a{x)dtu\'^ dtdx. 

for any solution of the damped equation ([3]). That means that there exists < C such that any 
solution of the damped equation satisfying ((75]) fulfills 

(76) E{u){T)<{l-C)E{um. 

Pick iV e N large enough such that {I - C)'^ C{Rq) < /2. 

Corollary 12.41 and ([74]) allow us to choose 6 small enough such that the assumption 

||(wo,^il)||£- < Rq] II('"0,'"i)|Il2x//-i < 5 

implies 

(77) \\{u{nT),dtu{nT))\\^,^H-^<e, < n < iV. 

So, with that choice, we have E{u){NT) < (1 — C)^ E{u)[{S). Then, by the energy decreasing, for 
any t > NT, we have 

\\{u,dtu)it)\\l,^j,,, < 2E{u){t) < 2E{u){NT) < e\ 
Therefore, the decay estimate (|76|) is true on each interval [nT, [n + l)r], n G N and we have 

E{u){nT) < (1 - CYE{u){Q) 
which yields the result. □ 
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Proof of Theorem W . 1\ Since the equation is reversible, we can assume {uq,ui) = (0,0). By a com- 
pactness argument, we can select a G Cq°(u;) with a{x) > r/ > for x inCo where Cj satisfies Assump- 
tion 10.21 We will first use the damping term a{x)'^dtu as a term of control. We apply Theorem 10.21 
and Theorem 13.21 once the energy is small enough. □ 

Acknowledgements. The author deeply thanks his adviser Patrick Gerard for drawing his attention 
to this problem and for helpful discussions and encouragements. 
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